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ABSTRACT 


The  work  presented  here  is  in  six  Sections.  In  the  first,  wc  review  the  main 
differences  between  the  piasma  environments  in  geostationary  orbit  and  low  polar 
orbit  with  regard  to  nigh- voltage  charging  situations.  In  Section  2,  we  develop  a 
simple  rough  estimate  of  the  required  conditions  for  overall  charging  of  a  large 
spacecraft  in  low-orbit  auroral-zone  conditions.  The  results  indicate  that  for  any 
given  spacecraft,  surface  potentials  are  likely  to  depend  more  strongly  on  the 
ratio  of  ambient  flux  of  high-energy  electrons  to  that  of  all  ions  than  on  any  other 
environmental  parameter,  and  this  prediction  has  been  corroborated  by  results  of 
Gussenhoven  et  al  for  the  DMSP  satellites.  In  Section  3,  we  present  results 
from  a  calculation  of  escape  currents  of  electrons  emitted  from  negatively- 
charged  spacecraft  surfaces  having  various  orientations  relative  to  the  direction 
of  the  local  magnetic  field  B.  The  suppression  of  such  currents  by  B  effects 
indicates  that  on  mostly-dieiectric  large  spacecraft  such  as  the  Shuttle  Orbitcr, 
local  charging,  especially  on  surfaces  nearly  parallel  to  B,  may  occur  in 
ionospheric  conditions  which  do  not  produce  overall  charging.  This  extends,  to 
lower  equivalent  temperatures,  the  range  of  auroral-electron  plasma  conditions 
in  which  one  can  expect  such  a  spacecraft  to  undergo  high-voltage  charging.  We 
aiso  investigate  surface  currents  of  nonescaping  emitted  electrons.  In  Section  4, 
we  propose  a  wake-induced- barrier-effect  mechanism  which  also  suppresses 
escape  of  emitted  electrons.  This  effect  appears  to  permit  beam-induced  as  well 
as  poiar-orbit  high-voitage  charging  to  occur.  In  Section  b,  we  present 
preliminary  results  of  numerical  simulation  work  directed  toward  making 
detailed  tests  of  the  predictions  made  in  Sections  2-4.  Section  6  contains  some 
concluding  remarks. 
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i.  INTRODUCTION  TO  THF.  AURORAL -ZONE  CHARGING  PROBLEM 


The  prediction  of  high-voltage  charging  or  other  environmental  effects  on  a 
spacecraft  in  low  Earth  orbit  is  more  difficult  than  in  geosynchronous  orbit,  for 
at  least  three  reasons:  (a)  space  charge  effects  (on  sheath  and  wake  potentials) 
are  more  important,  because  space-charge  densities  are  much  higher  (the  Debye 
length  is  no  longer  >>  typical  spacecraft  dimensions)  (b)  ion  flow  effects  are 
more  important,  because  spacecraft  orbital  speed  >  ion  thermal  speeds  (c)  the 
geomagnetic  field  B  is  likely  to  have  an  important  influence  on  charged-particle 
motions  because  B  is  now  much  larger,  and  not  all  of  the  average  particle 
gyroradii  of  importance  are  any  longer  >>  typical  spacecraft  dimensions.  In  the 
case  of  the  Shuttle  Orbiter,  othei  complications  arise.  The  most  important  one 
is  that  its  surface  is  almost  entirely  dielectric  rather  than  conductive.  The 
possibility  therefore  exists  that  local  charging  may  occur  on  the  Orbiter, 
especially  on  its  downstream-facing  surfaces,  in  ionospheric  conditions  which  do 
not  produce  overall  charging.  Another  complication  for  the  Orbiter  is  that 
charged -particle  mean  free  paths  may  not  always  be  much  larger  than  spacecraft 
size,  especially  during  thruster  firings  and  water  dumps.  A  further  complication 
arises  from  evidence  that  spontaneous  oscillations  occur  in  some  parts  of  the 
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Orbiter’s  ionized  wake.  Althouph  there  is  no  evidence  that  these  can 

O 

substantially  affect  the  likelihood  of  high  voit age  charging,  this  issue  is  largely 
unexplored.  Another  complicating  factor  is  the  geometrical  complexity  of  the 
Orbiter’s  cargo  bay  area.  When  this  study  was  begun,  polar  orbit  Shuttle  flights 
appeared  to  be  imminent,  and  its  primary  motivation  was  the  need  to  gain  an 
understanding  of  all  important  effects  governing  the  possibility  of  high-voltage 
charging  on  the  Orbiter’s  external  surfaces. 

Table  1  summarizes  some  important  characteristic  lengths  and  speeds  for 
low-orbit  conditions.  A  surprising  feature  of  this  Table  is  that  the  sheath 
thicknesses  indicated  are  much  larger  than  the  ambient  Debye  !  .ength,  but  this  is 
because  the  sheath  potentials  are  much  larger  than  the  ambient-particle  thermal 
energies.  These  distances  are  at  most  comparable  to  typical  spacecraft 
dimensions,  in  contrast  with  the  geosynchronous  situation.  Table  1  also  shows 
that  secondary  electrons  have  an  average  gyroradius  <<  typical  spacecraft 
dimensions,  so  their  escape  will  be  inhibited  strongly  on  surfaces  which  are 
.nearly  parallel  to  the  magnetic  field  B  (Fig.  3.1),  while  auroral  electrons  have 
an  average  gyroradius  )>  typical  spacecraft  dimensions,  so  their  collection  will 
be  affected  only  moderately,  except  for  very  large  spacecraft.  We  return  to  this 
question  in  Section  3.  Also  evident  from  Table  1  is  the  large  value  of  the  ion 
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speed  ratio  (spacecraft  speed/ion  most-probable  thermal  speed)  in  low-orbit 
conditions.  In  these  conditions,  ion  collection  on  downstream  surfaces  will  be 
inhibited.  If  a  surface  is  simultaneously  downstream  and  nearly  parallel  to  the 
magnetic  field,  as  is  likely  to  be  the  case  in  the  aurora:  zones,  then  the  tendency 
for  high-voltage  charging  to  occur  on  it  will  be  greatly  increased  (Fig.  3.2). 

High-voltage  charging  of  large  spacecraft  by  the  low-polar-orbit  plasma 
environment  was  first  predicted  by  Parks  and  Katz  (1981).  The  first  observation 
of  such  charging  was  made  by  Gussenhoven  et  al  (1985)  using  the  DMSP  6  and  7 
satellites.  Their  observations  indicate  that  high-voltage  charging  (conventionally 
defined  as  involving  spacecraft  surface  potentials  at  least  100  volts  different 
from  that  of  the  surrounding  plasma)  will  be  environmentally  produced  only  when 
a  spacecraft  encounters  an  "auroral"  electron  plasma  which  imposes  on  it  a 
sufficiently  large  ratio  of  "hot"  electron  ambient  flux  to  total  ion  (ambient  or 
ram)  flux  (see  their  Fig.  7  and  associated  discussion).  This  corroborates  a 
prediction,  made  by  Laframbois°  (1985b),  that  for  any  given  spacecraft,  surface 
potentials  are  likely  to  depend  more  strongly  on  the  ratio  of  ambient  flux  of  high- 
energy  electrons  to  that  of  all  ions,  than  on  eny  other  environmental  parameter; 
see  also  Laframboise  and  Parker  (1986,1987).  Section  2  of  this  Report 
contains  a  derivation  of  this  prediction.  The  theoretical  prediction  of  Parks  and 
Katz  (1981),  Katz  and  Parks  (1983),  indicates  that  such  charging  is  more  iikely 
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for  larger  spacecraft  because  electron  collection  increases  more  rapidly  with 
spacecraft  size  than  :on  collection  does. 

It  is  generally  accented  also  that  high-voltage  charging  will  occur'  only  when 
the  ambient  electron  plasma  has  a  relatively  large  fraction  of  its  total  flux  at 
energies  well  above  the  secondary-yield  maximum  of  the  spacecraft  surface 
material (s),  so  that  secondary  electrons  are  not  emitted  in  sufficient  numbers  to 
discharge  the  spacecraft.  This  situation  involves  "threshold"  behaviour,  with  a 
large  increase  in  surface  charging  produced  by  very  small  changes  in  the  ambient 
distribution  near  the  threshold  condition  (Rubin  et  ai,  1978;  Garrett  and  Rubin, 
1978;  Besse,  1981;  Garrett,  1981;  Laframboise  et  al,  1982;  Meyer-Vernet, 
1982;  Laframboise  and  Kamitsuma,  1983;  Lai  et  al,  1983;  Olsen,  1983; 

Mullen  et  al,  1986;  Katz  et  al,  (1986).  It  is  also  accepted  that  the  spacecraft 
must  be  in  darkness,  because  the  flux  of  photoelectrons  from  spacecraft  surfaces, 
which  is  in  the  range  10-50  uA/m2  for  most  surface  materials  (Feuerbacher  and 
Fitton,  1972;  Grard,  1973)  is  generally  larger  than  the  flux  of  high-energy 
auroral  electrons  encountered  during  known  polar-orbit  charging  events,  which  is 
usually  about  10  pA/m2  (Gussenhoven  et  ai,  1985,  Fig.  11).  However,  this  is 
not  always  the  case:  the  two  "tvpe-1"  auroral  -electron  distributions  reported  by 
Yeh  and  Gussenhoven  (1987,  Table  2a)  have  fluxes  of  57  and  42  pA/m2. 
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Laframboise  (1985a,  1988)  pointed  out  that  magnetic-field  effects  may 
suppress  secondary-electron  escape  and  thereby  increase  the  likelihood  of  high- 
voltage  charging.  Escape  of  photoelectrons  should  be  suppressed  similarly, 
because  secondary  electrons  and  photceiectroi  is  both  have  relatively  low  average 
energies  of  emission,  such  that  Maxwellian  fits  to  their  emission-velocity 
distributions  give  temperatures  of  about  3eV  and  1 .5eV  respectively  (Sternglass, 
1954;  Hachenberg  and  Brauer.  1959:  Chung  and  Everhart,  1974;  Hinteregger  et 
al,  1959;  Feuerbacher  and  Fitton,  1972;  Grard,  1973;  Wrenn  and  Heikkila, 
1973).  However,  no  examples  of  davlight  polar-orbit  high-voltage  charging  have 
so  far  been  seen.  In  Section  3  of  this  Report,  we  examine  this  suppression 
effect  in  detail.  In  Section  4,  we  examine  another  effect,  the  wake-induced- 
barrier  effect,  which  can  also  suppress  secondary-electron  or  photoelectron 
emission  from  surfaces  in  the  wake  of  a  large  dielectric-covered  spacecraft. 
Either  of  these  effects  may  produce  local  high-voltage  charging  in  circumstances 
which  would  net  produce  overall  charging,  in  particular,  they  extend,  to  lower 
equivale' '  .  •  nperatures,  the  range  of  auroral -electron  plasma  conditions  in 
which  one  car  >  xpect  such  a  spacecraft  to  undergo  high-voltage  charging. 
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In  Section  5  of  this  Report,  we  present  preliminary  results  from  a  numerical 
simulation  which  was  constructed  for  the  purpose  of  making  detailed  tests  of  the 
predictions  made  in  Sections  2-4.  This  work  involves  calculations  of  floating- 
potential  distributions  on  infinite  dielectric  cylinders  in  collisionless  plasma 
crossflows  whose  properties  model  those  of  the  auroral  plasma. 

Section  6  contains  some  concluding  remarks.  We  use  SI  units  throughout. 
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2.  ESTIMATE  OF  REOUIRED  ENVIRONMENTAL  CONDITIONS  FOR  LOW-POLAR- 
ORBIT  CHARGING. 

In  this  Section,  we  develop  a  simple  rough  estimate  of  the  conditions  necessary 
for  high-voltage  whole-body  charging  of  a  large  spacecraft  in  low  polar  orbit,  and  we 
show  that  spacecraft  surface  potentials  are  likely  to  depend  more  strongly  on  the 
ratio  of  ambient  flux  of  high-energy  electrons  to  that  of  oil  ions  than  on  any  other 
applicable  environmental  parameter.  To  do  this,  we  make  the  following 
approximations: 

(1)  In  this  Section,  we  assume  that  magnetic-field  effects  on  charged-particle 
motion  are  negligible.  This  assumption  should  be  acceptable  for  initial  estimates 
because  the  gyroradii  of  ions  and  high-energy  electrons  are  generally  a  few  meters 
or  larger,  especially  in  a  high-voltage  sheath  (Table  1),  and  collection  of  "cold" 

(-'-0.1  eV)  ionospheric  electrons  by  a  negatively-charged  spacecraft  will  be  very 
small,  so  their  density  is  well-approximated  by  a  Boltzmann  factor,  independently  of 
the  presence  of  a  magnetic  field.  In  Sec.  3,  we  return  to  the  question  of  magnetic- 
field  effects. 

(2)  We  assume  that  ambient  high-energy  electrons  have  an  isotropic  velocity 
distribution.  Large  departures  from  this  have  been  observed  in  auroral-plasma 
conditions  (Fennell  et  al,  1981;  Lin  and  Hoffman,  1982;  W.J.  Burke,  1984,  private 
communication;  Yeh  and  Gussenhoven,  1987),  but  this  should  not  seriously  affect  the 
type  of  rough  estimate  made  here.  Parks  and  Katz  (1981)  and  Katz  and  Parks 
(1983)  assumed  both  the  ion  and  electron  fluxes  to  be  unidirectional;  we  discuss  this 
point  later  in  this  Section. 
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(3)  We  ignore  secondary-electron  emission:  magnetic-field  effects  would  tend  to 
suppress  this  on  some  parts  of  the  spacecraft  in  any  case  [Laframboise  (1983a, 
1985);  Sec.  3]. 

(4)  We  assume  that  the  spacecraft  is  a  unipotential  sphere,  large  compared  to  the 
typical  ambient  Debye  length  of  ^  1  cm  (Table  1).  We  consider  only  overall  charging 
of  the  spacecraft.  This  neglects  the  possibility  that  local  high-voltage  charging  may 
occur,  especially  on  surfaces  in  the  spacecraft  wake  (Sections  3-5). 

(5)  We  assume  that  both  ions  and  electrons  have  double-Maxwellian  velocity  distri¬ 
butions,  with  the  colder  component  in  either  case  having  a  temperature  of  0.1  eV, 
and  the  hotter  1  keV  or  larger.  In  the  spacecraft  reference  frame,  these  are  super¬ 
posed  on  a  drift  velocity  equal  and  opposite  to  the  spacecraft  velocity. 

(6)  Ions  are  assumed  to  be  either  H+  or  0+  . 

Note  that  assumption  (3)  could  cause  a  false  prediction  that  high-voltage  charging 
occurs,  while  assumption  (4)  could  cause  a  false  prediction  that  it  does  not.  The 
effects  of  assumptions  (1),  (2),  and  (5)  are  less  clear;  these  could  conceivably  either 
increase  or  decrease  predicted  surface  potentials.  With  regard  to  (6),  assuming  that 
the  ions  are  H*  results  in  maximum  wake-filling  by  ions.  If  there  are  any  electri¬ 
cally-isolated  surfaces  in  the  spacecraft  wake,  this  would  result  in  decreased  surface 
potentials  (magnitudes);  assuming  0+  gives  the  reverse  (see  also  Sec.  4). 

Probably  the  most  serious  difficulty  in  formulating  a  theory  for  low-orbit 
charging  is  the  prediction  of  ion  collection  on  downstream  surfaces.  As  mentioned  in 
assumption  (4)  above,  we  avoid  this  difficulty  by  considering  only  total,  rather  than 


8 


local,  ion  collection,  on  a  unipoter.tial  sphere.  We  return  to  the  question  of  local  ion 
collection  in  Sec.  4.  Kanal  [1962,  F.q.  (63)]  gives  an  expression  for  the  total  ion 
current  collected  by  such  a  sphere  from  a  drifting  Maxwellian  plasma  in  the  limit  of 
zero  potentials  (relative  to  space  potential),  as  follows: 


ij  =  £  M  (Sj  + 


'  er'f (S . }  +  exp(  S  )] 

.  TD  .  1  1 


(2.1) 


where  i.  =  I./I  I  .is  the  ion  random  current  en.  (kT./2rrm.)  S.  = 
l  r  oi  oi  loo  i  l  i 

U/(2kT./m.)^  is  the  ion  speed  ratio,  U  is  the  ion  drift  speed  relative  to  the  space¬ 
craft,  e  is  the  magnitude  of  the  electronic  charge,  k  is  Boltzmann’s  constant,  and  rrn, 
T.,  and  n.^  are  ion  mass,  temperature,  and  ambient  number  density.  We  assume 
that  U  =  8  km/sec,  corresponding  to  low  circular  orbit. 


We  need  to  take  account  of  the  effect  of  a  large  ion-attracting  surface  potential  on 

ion  collection,  in  the  limit  of  small  Debye  length  A^  compared  to  the  sphere  radius 

r  .  To  do  this,  we  use  a  result  of  Parrot  et  al  (1982).  These  authors  show  that  for 
s 

a  probe  in  a  collisionless,  nonmagnet izod,  Maxwellian  plasma  having  T./Tg  =  1  and 
without  ion  drift,  and  in  the  limb  when  An/r  — >  0  but  -e<p  /kT  >>  1  [where  <p  is 
surface  potential  relative  to  space,  and  these  limits  must  be  approached  in  such  a 
way  that  (-e$s/kT)  (A^/r^4'3  remains  <<  1,  i.e.,  sheath  thickness  remains  <<  sphere 
radius],  the  ion  (att.racted-particle)  current  is  larger  than  the  random  current  by  a 
factor  of  1.45.  This  factor  represents  the  effect  of  "presheath"  electric  fields  on  ion 
collection.  Even  though  several  of  their  assumptions  are  unfulfilled  in  our  case,  the 
resulting  effects  on  ion  collection  are  probably  small  enough  for  our  purposes.  We 
therefore  multiply  F.q.  (2.1)  by  the  same  factor  to  obtain  an  estimate  of  total  ion 
collection  as  Influenced  by  surface-potential  effects.  The  resulting  ion-current 
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dependence  on  ion  speed  ratio  is  plotted  in  Fig.  2. 1  .  For  O*  ions  at  T  =  0.1  eV 
(.1  160K),  H+  at  0.1  eV,  O'*  at  IkeV,  and  IT  at  IkeV,  we  have  S.  =  7.31,  1.83, 
Q.0731,  and  0.0183  (the  latter  two  are  effectively  zero),  respectively.  The  corres¬ 
ponding  ion-current  enhancement  factors  (values  of  i.)  from  Fig.  2.1  are  9.50,  2.69, 
1.45,  and  1.45,  respectively. 

If  the  ambient  ions  are  rT,  the  ion  collected  current  is  now  given  by: 


I.  =  47m  2 
1  s 


en 


r 

ic{  2 


k  i  . 


_ic 

7rm. 

i 


(2.69) 


47m 


en. 


ih 


P.Trrm 


\  (1-45) 


(2.2) 


where  the  subscripts  ic  and  ih  refer  to  the  cold  and  hot  ion  populations.  If  the  ions 
are  0+,  then  the  factor  2.69  in  (2.2)  should  bo  replaced  by  9.50. 

The  electron  collected  current  is: 


If  high-voltage  charging  occurs,  then  -ecp^  >>  kT^,  and  the  first  term  on  the  right- 
hand  side  of  tins  equation  becomes  negligible. 
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For  current  balance,  I.  =  I  .  This  leads  to: 

’  :  e 

2'69"ic^TT  +  1-‘15nlhV'TT  “  neh  Vm/Sv  'TT^e  "e,'*’s!/tcTeh 


(2.4) 


where  V^m./m  =  43  for  hT  ions.  Therefore: 
v  e 


el<tsIATeh=  in 


43nehVTh 


2.69n.  VT.  +  1.45n.,  VT., 
ic  -  ic  lh  -  lh 


i?..5) 


for  H*  ions,  with  43  and  2.69  replaced  by  172  and  9.50  for  O*  ions.  This  is 
equivalent  to: 


e|*sl/kTeh  "  ,n 


hot-electron  ambient  flux 


2.69  (cold-ion  ambient  flux)  +  1.45  (hot-ion  ambient  flux) 


(2.6) 

For  high-voltage  charging  to  become  probable,  the  argument  of  the  in  function 
must  be  close  to  or  larger  than  e  2  2.72,  :.e: 


hot-clcctron  ambient  flux 


2.69 (coid-ion  ambient  flux)  +  1 .45 (hot-ion  ambient  flux) 


£  2.72. 


(2.7) 


For  0+/H+  mixtures  and  for  hot-ion  temperatures  other  than  1  keV,  generalization 
of  this  result,  is  straightforward.  Since  any  hot  ions  are  likely  to  have  WTic~  10‘- 
the  hot-ion  ambient  flux  will  exceed  the  cold-ion  ambient  flux  if  the  hot  ions 
constitute  more  than  about  1%  of  the  total  ambient-ion  number  density,  equation 
(2.7)  indicates  that  the  onset  of  high-voltage  charging  can  be  expected  to  depend 
primarily  on  the  ratio  of  hot-electron  ambient  flux  to  the  ambient,  flux  of  all  ions,  as 
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mentioned  at  the  beginning  of  this  Section.  This  completes  our  argument  in  support 
of  this  conclusion. 

In  analyzing  spacecraft  data,  one  is  therefore  likely  to  find  better  correlation  of 
spacecraft  voltages  with  the  ratio  which  appears  in  F.q.  (2.7),  or  something  nearly 
equal  to  it,  than  with  any  other  measurable  quantity,  such  as  electron  or  ion  density 
or  average  energy,  taken  individually.  This  expectation  has  been  borne  out  in  recent 
work  by  Gussennoven  et  al  (1985,  Fig.  7),  involving  charging  data  from  the  DMSP 
F6  and  F7  satellites.  In  calculating  values  of  this  ratio,  the  ambient  fluxes  which 
are  involved  need  to  have  been  measured  simultaneously  on  the  same  spacecraft. 

Tver,  though  the  approximations  made  in  deriving  (2.7)  are  severe,  and  the  precise 
dependence  of  spacecraft  voltages  on  this  ratio  may  therefore  differ  substantially 
from  that  given  in  Eq.  (2.7)  (and  the  coefficients  in  (2.7)  will  need  to  be  modified  if 
FT  dominates),  our  general  conclusion,  i.e.  that  spacecraft  voltages  should  correlate 
most  strongly  with  this  ratio,  or  something  nearly  equal  to  it,  is  likely  to  remain 
valid.  Furthermore,  the  dependence  of  spacecraft  voltages  on  this  flux  ratio  is  likely 
to  retain  an  approximately  exponential  form.  In  situations  where  most  secondary  and 
backscattered  electrons  emitted  by  the  spacecraft  will  escape  (see  Sec.  3),  primary- 
electron  incident  fluxes  will  be  approximately  cancelled  for  many  spacecraft 
materials  by  electron  escape  at  incident  energies  up  to  a  few  keV  (Laframboise  et 
si,  1982a,b;  Laframboise  and  Kamitsuma,  1983;  Lai  et  al,  1983),  so  the  hot- 
electron  ambient  flux  term  in  (2.7)  needs  to  be  modified  accordingly. 

The  most  serious  approximation  made  in  deriving  (2.7)  is  probably  item  (4)  in  the 
list  at  the  beginning  of  this  Section.  This  is  because  ion  fluxes  on  downstream 
surfaces  are  iikeiy  to  be  very  much  smaller  than  their  average  over  the  entire 
spacecraft  (see  also  Sections  3-5).  They  are  also  likely  to  be  strongly  dependent  on 
spacecraft,  geometry,  local  surface  potential  distribution,  and  0+/H+  concentration 
ratio.  Therefore,  the  critical  value  of  ambient  flux  ratio,  at  which  the  onset  of  high- 
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voltage  charging  occurs,  is  likely  to  vary  substantially  among  spacecraft  having 
different  geometries  and  surface  materials.  In  particular,  for'  spacecraft  having 
electrically-isolated  downstream  surfaces,  this  critical  ratio  is  likely,  because  of 
local  charging  on  these  surfaces,  to  be  much  lower'  than  for  spacecraft  which  have  an 
entirely  conductive  surface  (Sections  3  5), 

Furthermore,  in  contrast  with  the  situation  for  total  ion  collection,  no 
known,  simple,  reliable  method  has  been  available  for  estimating  ion  fluxes  on 
downstream  surfaces.  Parks  and  Katz  (1983,  Sec.  4;  1985)  have  developed  an  ior, 
flux  and  density  calculation  for  the  downstream  point  on  a  sphere  in  a  model  potential 
which  has  a  given,  simple  analytic  form.  In  Sen.  4.3  and  Appendix  C  we  use  a 
similar  method  to  develop  an  ion  density  calculation  for  trie  downstream  location  or.  a 
cylinder  in  a  collisionless  crossflow.  Detailed  numerical  simulation,  which 
includes  realistic  self-consistent  spacecraft  sheath  potential  distributions,  and  which 
probably  needs  to  involve  at  least  some  ion  orbit-following,  will  be  essential  to 
determine  the  accuracy  of  such  approximations.  In  Sec.  5,  we  report  on  preliminary 
results  from  a  calculation  of  this  kind. 

So  far,  we  have  not  mentioned  the  difficulties  which  can  arise  in  measuring  the 
ambient  ion  fluxes  which  appear  in  Eq.  (2.7).  So  far,  we  have  also  defined  "ambient 
flux"  to  be  that  measured  in  an  Earth-fixed  reference  frame.  The  alternative  would 
be  to  define  it  as  that  measured  in  the  spacecraft  frame,  i.e.,  including  ram  effects. 
Ion  fluxes  measured  by  spacecraft  instruments  are  strongly  influenced  by  ram 
effects.  In  fact,  the  numerical  factors  2.69,  1.45,  and  9.50,  which  appear  in  Eq. 
(2.7)  and  the  associated  discussion,  already  constitute  a  rough  ram-effect  correction, 
but  for  total  current  to  a  sphere,  not  for  local  collection  by  a  forward-facing  instru¬ 
ment  aperture.  It  may  happen  that  the  ram-effect  correction  factors  for  an 

instrument  are  nearly  equal  to  the  above  factors,  so  that  trie  instrument  measure- 
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ment,  without  any  correction,  already  gives  a  good  estimate-  of  the  denominator  of  F.q. 
[2.7).  In  any  case,  the  response  of  the  instr  ument  will  depend  on  its  geometry,  and 


this  problem  has  already  been  treated  by  other  authors  (Parker,  1970;  Parker  and 
Whipple,  1  970;  Whipple  et  al,  1974;  Chang  et  al,  1979;  Singh  and  Baugher,  1981; 
Comfort  et  al,  1987;  l.aframboise,  1983b),  so  we  do  not  discuss  it  here. 

Parks  and  Katz  (1981)  and  Katz  arid  Parks  (1983)  have;  estimated  charging 

poterriais  or.  snherU.cil  spacecraft  of  0.5m  and  3m  radius,  assuming  that  the  ions  are 

O*,  th.e  hot-electron  temperature  T  ^  is  5  keV,  and  spacecraft  speed  is  8  km/scc. 

Their  results  can  be  compared  directly  wi*h  those  given  by  our  F.qs.  (2.5)  -  (2.7). 

They  have  used  the  thecry  of  Langmuir  and  Blodgett.  (1924)  to  obtain  values  for 

sheath  radius  as  a  function  of  spacecraft  potential.  They  present  spacecraft 

potentials  as  functions  of  the  ratio  k  of  hot  ("precipitating")  electron  ram  current  to 

ion  ram  current.  To  make  a  comparison,  their  value  of  k  needs  to  be  expressed  in 

terms  of  our  ambient  flux  ratio.  They  have  assumed  the  ambient  electron  flux  to  be 

unidirectional.  To  convert  to  an  equivalent  isotropic  flux,  we  note  that  current  to  a 

sphere  =  Anr  2  x  isotropic  (random)  flux,  but  =  nr  2  x  undirectional  (ram)  flux, 
s  s 

Therefore,  equivalent  isotropic  flux  =  {  x  unidirectional  flux,  for  a  sphere. 


Also  for  a  sphere,  the  ratio  of  ion  ram  to  random  currents  is 

U/(8kT./rrm.)^  =  \  VTrS..  Using  S.  =  7.31,  this  ratio  ~  6.48,  so  therefore: 
l  i  -  i  i 

y  _  hot  electron  ram  current _ 

6.  48  x  total  ion  random  current 


nr  2  x  I  lot  electron  ram  flux 
s _  _ _ 

6 .48  x  Anr  7  x  total  ion  random  flux 

s 

1  hot  electron  (equivalent)  random  flux 

6.48  total  ion  random  flux 


(2.8) 


x  our  flux  ratio  R. 


With  coefficients  for  0+  used,  our  F.q.  (2.6)  gives: 
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<*>  =  -5000  Jn  (R/9.50).  (2.9; 

Figure  2.2  shows  our'  resuh  and  theirs  [from  their  Fig.  3  (1981)  or  Fig.  2  (  19S3)j, 
plotted  together.  At  larger  potentials,  the  combined  set  of  results  shows  a  monotonic 
progression  toward  increased  charging  for-  larger  spacecraft.  For  <  350V,  their 
5rr  sphere  shows  more  charging  than  our  large-radius- limit  sphere.  This  is  because 
•heir  ion -current  enhancement  factor,  which  is  determined  hv  the  size  of  a  sharp- 
edged  Langmuir-Blodgett  sheath,  falls  below  ours,  which  includes  the  effect  of  a 
quasineutral  presheath.  This  discussion  suggests  that  the  tendency  toward  high- 
voltage  charging  always  increases  with  spacecraft  size,  but  magnetic-field  effects 
may  modify  this  ([  .aframboise,  1983a,  Sec.  1;  !  .aframboise,  1988,  Sec.  1;  see  also 
Sec.  3.1  of  this  Report).  The  corresponding  curves  for  local  charging,  on  surfaces 
in  a  spacecraft  wake,  will  lie  to  the  left  of  those  shown  in  Fig.  2.2,  but  these 
remain  to  be  computed  numerically,  as  we  have  done  for  a  particular  spacecraft 
geometry  in  Sec.  5. 
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CALCULATION  OF  ESCAPE  CURRENTS  OF  ELECTRONS  EMITTED  FROM 
NEGATIVELY -CHARGED  SPACECRAFT  SURFACES  IN  A  MAGNETIC 
FIELD. 

TNTRO.:  ■!  'CTION  TO  THE  ESCAPE-  'URRENT  PROBLEM 

Ir.  Sec.  1 ,  it  was  pointed  out  that  ir.  lo w  Earth  orbit,  the  geomagnetic  field 
B  ir  core  likely  tc  Evvt  an  important  ;•  fh.ence  on  charged -particle  mo* ions 
than  ir  goosynchronc  us  orbit,  because  B  is  now  much  stronger,  and  not  all  of 
‘he  average  particle  gyroraaii  of  importance  are  any  longer  >>  typical 
spacecraft  dimensions.  In  this  Section,  we  wish  to  investigate  an  important 
consequence  of  this  situation,  which  concerns  the  escape  of  electrons  emitted 
from  spacecraft  surfaces.  The  emitted  velocity  distributions  of  secondary 
electrons  and  photoelectrons  are  frequently  approximated,  in  both  speed  and 
direction,  by  Maxwell lans  corresponding  to  a  temperature  T  close  to  3  eV  for 
secondary  electrons  (Stemgiass,  1  959;  Harbenberg  and  Brauer,  1959;  Chung 
and  Everhart,  1979)  and  1.5  eV  for  photoelectrons  (Hir.teregger  et  al,  1959; 
Eeuerbacher  and  Fitton,  1972;  Grard,  1973;  Wronr,  and  Heikkila,  1973). 
Backscattered  electrons  have  energy  distributions  which,  generally  bear  less 
resemblance  to  Max  we!  liars,  and  they  also  have  much  larger  average  energies 
(Stemgiass,  1959,  Fig.  8),  so  the  treatment  to  be  developed  herein  cannot  be 
applied  as  readily  to  them.  However,  their  contribution  to  the  prevention  of 
high-voltage  spacecraft  charging  is  generally  less  decisive  than  that  of  either 
secondary  electrons  or  photoelectrons  because  their  current  contribution  by  its 
nature  is  always  smaller  than  that  of  the  incident  "primary"  electrons.  A 
review  of  the  properties  of  electron  emission  from  surfaces  in  space  has  been 
given  by  Whipple  (1981). 
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In  low  Earth  orbit,  in  the  auroral -zone  geomagnetic  field  (!B!=0.44  gauss 
=  4.4  x  10'5T),  the  gyroradn  of  a  "typical"  3eV  secondary  electron,  a  1.5 
eV  photoeieetron,  and  a  10  keV  auroral  electron  are  13  cm,  9  cm,  and  8 
m,  respectively.  The  average  gyroradius  of  "cold"  ionospheric  electrons 
(temperature  .  .1.1  eV)  m  the  same  B  is  even  smaller  (2  cm),  but  this  is 

not  an  important  parameter  in  most  cases  because  these  electrons  are  receded 
if  the  spacecraft  potential  is  negative,  and  their  density  is  then  well- 
approximated  by  a  Boltzmann  factor,  which  is  unaltered  by  B  effects. 

The  reason  whv  B  affects  escape  of  emitted  electrons  is  shown  in  Fig.  3.1. 
In  Fig.  3. Hat ,  the  spacecraft  surface  is  perpendicular  to  B,  and  the  emitted 
electrons,  which  experience  an  electric  force  -eE  directed  away  from  the 
surface,  all  escape,  helping  to  discharge  it.  In  Fig.  3.1(b),  the  spacecraft 
surface  is  nearly  parallel  to  B,  and  almost  all  of  the  emitted  electrons  return 
to  it,  even  ‘hough  they  stili  experience  an  electric  force  directed  away  from  it. 

These  electrons  ‘herefore  are  unable  to  help  discharge  it,  so  n  surface  nearly 
parallel  to  B  is  more  likely  to  charge  *o  a  large  negative  voltage.  Note  that 
‘he  component  of  E  which  i  i  perpendicular  to  B  results  only  in  an  E  x  B  drift 
oarallei  to  the  surface. 


For  any  object  much  large:-  than  13  cm,  the  escape  of  secondary  electrons 
and  photoelectrons  will  be  strongly  affected  by  this  process.  For  example, 
most  surfaces  on  the  Shuttle  are  effectively  "infinite  planes"  by  this  criterion. 
Or  the  other  hand,  the  average  gyroradius  of  high-energy  auroral  electrons  is 
comparable  to  Shuttle  dimensions,  so  the  deposition  of  these  e’ectrons  onto 
Shuttle  surfaces  is  likely  to  be  only  moderately  inhibited. 
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For  a  larger  object  (sire  >)■  8  r.  deposition  of  auroral  electrons  will 
also  become  strongly  orientation-dependent,  with  both  collection  and  escape  of 
electrons  now  being  inhibited  on  surfaces  nearly  parallel  to  B.  This  suggests 
that  high-voltage  charging  of  such  sur  faces  may  be  mere  likely  on  objects  of 
intermediate  size  than  on  either  larger  or  smaller  ones.  In  the  calculation  of 
Parks  and  Katz  11981),  Katz  and  Parks  (1983),  the  tendency  toward 
high-voltage  charging  increased  with  spacecraft  size  because  in  their  model, 
ion  collection  increased  less  rapidly  with  spacecraft  size  than  did  electron 
collection.  To  determine  which  of  these  two  effects  predominates  will  require 
more  detailed  calculations  than  have  been,  done  so  far.  The  theory  of  magnetic- 
field  effects  on  current  collection  by  objects  in  space  plasmas  is  still  very 
incomplete  (Parker  and  Murphy,  1967;  Linson,  1969;  Sanmartin,  1970; 
Laframbolse  and  Rubinstein,  1976;  Rubinstein  and  La  framboise,  1978,  1982, 
1983;  McCov  et  al,  19801. 
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As  already  mentioned,  strong  ion  flow  effects  also  are  generally  present  in 
low  orbit;  the  ion  speed  ratios  (flow  speed/most  probable  ion  thermal  speed) 
for  H  at  1  keV,  H1"  at  0.1  e'v’,  and  G  at  0.1  eV  are  0.02,  1.8,  and  7.3, 
respectively.  Whenever  the  latter  is  the  predominant  ion  species,  ion 
collection  on  downstream  surfaces  will  be  strongly  inhibited.  If  a  surface  is 
simultaneously  downstream  and  nearly  parallel  to  B,  as  is  likely  to  be  the 
case  m  the  auroral  zones,  then  ‘he  tendency  for  high-voltage  charging  to  occur 
on  will  be  greatly  increased  (Fig.  3.2).  An  extensive  literature  exists  on 
wakes  of  saielliVs  z  low  Perth  orbit.  Reviews  of  it  have  been  given  by 
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'*.986),  and  Stone  and  Samir  '1986).  Mont  of  this  literature  treats  only  low- 
voltage  wakes  and  nonduet ive  spacecraft.  Exceptions  include  Parks  and  Katz 
11  983a, bK  Parker  (1983)  and  Lafrarnboise  and  Parker  (1987). 

To  'straightforwardly"  include  B  effects  on  electron  escape  in  a  large  two 
or  three  dimensional  simulation  program  would  involve  the  numerical  integra¬ 
tion  of  very  large  numbers  of  emitted-eleetron  orbits.  The  resuiting  computing 
costs  usually  would  be  formidable,  especially  since  most  of  these  orbits  would 
have  relatively  large  curvatures.  A  desirable  alternative  is  to  "parameterize" 
the  situation  by  treating  in  advance  a  simplified  but  still  sufficiently  realistic 
model  problem.  Such  a  treatment  also  serves  several  other  purposes.  For 
particular  situations,  it  oermits  simole  estimates  of  the  effects  considered.  It 

V  'll 

may  help  to  stimulate  experimental  investigations  of  these  effects.  It  also 
permits  gaining  an  understanding  of  some  of  the  major  effects  which  govern 
the  overall  problem,  and  this  can  be  its  most  imDortant  use.  In  Sections 

i  •  i 

3. 2-3. 5,  we  develop  a  "parameterization"  of  magnetic-field  effects  on  esoapir.g- 
electror;  currents.  In  Sec.  3.6,  we  propose  a  simple  aDproximate  method  for 
calculating  the  space-charge-density  distribution  of  escaping  electrons. 
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3.2.  THEORY  FOR  ELECTRIC  FIE!!)  NORMAL.  TO  SURFACE 


in  order  to  develop  a  simplified  tr  eatment  of  electron  escape,  we  make  the 
following  approximations.  (1)  We  assume  (Sec.3.1)  that  electrons  are 
emitted  from  tire  spacecraft  surface  with  a  Maxwellian  velocity  distribution 
corresponding  to  a  temperature  T.  (2)  We  assume  that  the  spacecraft  surface 
is  an  infinite  plane,  and  the  electric  and  magnetic  fields  FI  and  B  outside  it 
are  uniform  and  time- independent.  (3)  In  this  Section,  we  assume  that  the 
electric  force  -eE  on  electrons  is  directed  along  the  outward  normal  to  the 
surface;  here  e  is  the  magnitude  of  the  elementary  charge.  This  assumption 
is  relaxed  in  Sec.  3.4,  in  order  to  permit  variations  of  potential  along  the 
surface  to  be  taken  into  account. 

If  -eE  is  normal  to  the  surface,  our  assumptions  imply  that  the  ratio  i  = 
I/I  of  escaping  to  emitted  flux  is  a  function  of  two  parameters:  the 
angle  6  between  the  surface  normal  and  the  direction  of  B  (Fig.  3.3),  and  a 
parameter  describing  the  strength  of  E.  A  convenient  choice  for  this  para¬ 
meter  is  the  difference  in  potential  across  the  average  emitted-electron 
gyroradius  a  =  (1/eB)  (irmkT/2)^, divided  by  kT/e,  where  m  is  electron  mass 
and  k  is  Boltzmann’s  constant. 

This  quotient  is: 

e  -Tf  (3,1) 

where  F.  a  |E|  and  B  =  j B  i . 
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This  quantity  also  has  an  alternative,  more  useful  interpretation:  it  is  the 
ratio  of  the  magnitude  |E  x  Bl./B3  of  the  F.  xB  drift  speed,  to  one-half  the 
average  thermal  speed  (8kT/7r  m)‘  of  the  emitted  electrons.  It  is  useful  to 
estimate  the  value  of  e  for  a  high-voltage  spacecraft  sheath  in  low-orbit 
conditions.  To  do  this,  we  use  the  spherical ly-symmetric  sheath  solution  of 
Al’pert  et  al  (1965,  Table  XXIV  and  Fig.  72).  For  a  sphere  of  radius  3m 
having  a  surface  potential  of  - 1  kV  or  -5  kV  relative  to  space,  in  a  collision¬ 
less  plasma  having  an  ambient  ion  temperature  of  0.1  eV,  number  density  of 
3  x  10s  cm-3,  and  resultant  (ion)  Debye  length  of  0.43  cm,  their  results  give, 
respectively,  sheath  thicknesses  of  2.6  and  6.1  m,  and  surface  electric  fields 
E  =  0.86  and  2.9  kV/m.  Using  B  =  4.4  x  10 ~ST  and  T  =  3  eV  (1.5  eV)  for 
secondary  electrons  (photoelectrons)  ,  we  then  obtain  e  =  33.9  and  1  14.2 
(47.9  and  161.5).  These  are  relatively  large  values,  whose  significance  can 
be  understood  if  we  consider  what  would  happen  if  e  were  infinite. 

In  this  limit,  it  is  easy  to  show  that  emitted  electrons  would  all  escape 
unless  B  were  exactly  parallel  to  the  surface  (0  were  90°)  .  This  can  be 
shown  as  follows.  In  this  limit,  emitted  electrons  would  have  no  "thermal" 
motion.  The  (y,z)  projection  of  their  motion  would  then  be  similar  to  that 
shown  in  Fig.  3.4.  This  motion  would  be  the  sum  of:  (i)  an  E  x  B  drift  in  the 
y  direction  (ii)  a  uniform  acceleration  along  B  ,  whose  projection  in  the  (y,z) 
plane  would  be  upward  (iii)  just  enough  gyromotion  to  produce  a  cycloidal  path 
when  combined  with  (i),  so  that  in  the  absence  of  (ii),  the  electron  would  (just) 
return  to  the  surface  at  the  end  of  each  gyroperiod.  In  the  presence  of  (ii), 
these  ''return  points"  are  displaced  upward  by  progressively  increasing  amounts 
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(Fig.  3.4),  so  the  electron  can  never  return  to  the  surface,  unless  B  is  exactly 
parallel  to  the  surface,  so  that  the  upward  component  of  -eE  along  B  vanishes. 

This  result  suggests  that  for  large  finite  values  of  e  (including  the  values 
calculated  above),  electron  escape  is  likely  to  be  almost  complete  except  for  0 
very  near  90°,  where  it  should  drop  to  zero  very  steeply.  The  occurrence  of 
high-voltage  charging  in  marginal  circumstances  may  therefore  depend  very 
strongly  on  the  precise  orientation  of  a  surface. 


The  escaping  electron  flux  is  given  by: 
'  =  ///  f(v)  H(v)  vo2  d3V0 
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where:  is  the  initial  velocity  of  an  emitted  electron,  f(vQ)  =  d3n/d3vQ  is 

the  velocity  distribution  of  emitted  electrons,  n  is  a  reference  number  density, 
and  H(v  )  is  equal  to  1  for  escaping  electrons  and  0  for  those  which  return  to 
the  surface.  The  emitted  flux  is: 

IQ  =  n(kT/27rm)*.  (3.3) 

We  also  introduce  the  dimensionless  velocity: 

u  =  v  (m/2kT)^  .  (3.4) 
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Equation  (3.2)  then  becomes: 
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which  is  in  a  form  suitable  for  numerical  summation.  The  quantities  u^.^  ^ , 

u, .  u  ,  are  the  values  of  u  for  which  H  changes  between  0 

hm,2  ’  hm,k _  oz  b 

max 

and  i  for  each  u  and  u  .  These  values  must  be  found  by  numerically 
ox  oy 

determining  which  particle  orbits  reimpact  the  surface.  These  orbits  can, 
however,  be  determined  in  analytic  form,  with  time  as  a  parameter.  To  do 
this,  we  use  the  coordinate  system  shown  in  Fig.  3.3,  together  with  a  y-axis 
(not  shown)  directed  into  the  plane  of  the  Figure.  The  equation  of  motion  for 
an  electron  is: 


V  =  -^(E+vxB). 


We  solve  this  with  the  initial  conditions  ^  =  y  =  n  =  0,v^ 
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(3.7) 


^  'V 

x  =  x/a,  y  =  y/a,  etc; 
r  =  wct  =  (eB/m)t. 

In  the  present  discussion,  r^and  ar'e  both  zero,  but  for  later  use  (Sections 
3.4  and  3.5),  we  have  retained  these  quantities  in  the  formulas  below.  We 
obtain: 

u  .  =  u  sin  0  +  'j  cos  0; 
o£  ox  oz 
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Equations  (3.8)  can  also  be  differentiated  to  find  dz/dr.  The  numerical 

procedure  for  finding  the  quantities  u^ljT  ^  in  Eq.  (3.5)  then  involves 

calculating  z  and  dz/dr  at  a  succession  of  points  along  an  orbit  (the  electron 

will  reimpact  during  the  first  gyroperiod  0  <  r  ^  2tt  if  at  all,  so  this  interval 

always  suffices),  and  making  the  appropriate  tests  on  these  quantities  to  find 

out  whether  the  orbit  reimpacts  or  escapes.  For  each  u  ,  and  u  this  is 

done  for  a  succession  of  values  of  u  .  These  tests  also  yield  the  local 

oz  1 

minimum  of  z(r)  if  one  exists.  Whenever  a  change  occurs  between  no  escape 
and  escape  from  one  such  value  of  uq2  to  the  next,  an  interpolation  using  these 
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minima  can  be  used  to  provide  the  corresponding  value  of  u, .  »  .  In  cases 

1 1  IT  1  5  lx 

where  thev  are  unavailable,  the  arithmetic  mean  of  the  two  successive  u 
values  is  used. 

We  have  chosen  the  abscissas  u  .  and  u  .  in  Eq.  (3.5)  by  first  solving 

ox,i  oy,j  1/7  b 

numerically  the  transcendental  equation  ^  +  £  erf(u^/N/2.)  =  k/n  for 

k  =  1,2, _ ,  n- 1 .  The  resulting  values  u^  then  subdivide  the  normalized 

Gaussian  distribution  exp(-u2)/Vrr  into  n  "slices"  whose  areas  decrease  as 
expl-^u^2)  when  n  is  large.  In  terms  of  these  u^,  we  then  choose  values 
uq  k  at  the  centroids  of  these  slices,  and  weights  c^  equal  to  the  areas  under 
each.  We  obtain: 


exp(-uk  j2)  -  exp(-uk2) 

u  - - 

°‘  Vrr(erf  u^  -  erf  u^ 

ck  =  i(erf  uk  -  erf  u^) 


(3.9) 


(3.10) 


for  k  =  1,2, _ ,n.  A  convenient  method  for  calculating  the  required  values  of 

erf  u  =  1-erfc  u  has  been  given  by  Shepherd  and  Laframboise  (1981).  The 
resulting  values  u  ,  are  then  used  to  provide  the  required  values  of  u  . 
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and  u  ,  and  the  c.  are  used  to  provide  values  of  ( 1  /Vrr)  Au  .  exp(-u  . 2) 
oy,j  k  ox  ^  1  ox^i 


and  (l/V7r)Au  .  exp(-i 


i,  for  use  in  Eq.  (3.5).  We  have  provided  values 


of  uQZ  for  use  in  determining  the  u^m  k  values  in  Eq.  (3.5)  by  solving  the 
equation  l-exp(-u  ,2/2)  =  k/n  for  k  =  0,1,2,...,  n-1.  This  gives  u  ,  ~ 

0Z,K  OZ^K 

(  2in[  1/(1 -k/n)]}^;  these  values  are  distributed  most  densely  near  uqz  =  0,  but 
still  densely  enough  at  large  uqz  that  the  resulting  intervals  give  vanishing  flux 

contributions  in  this  limit.  This  completes  the  definition  of  the  procedure  used 
for  calculating  the  ratio  I/I  of  escaping  to  emitted  flux.  The  computer 
program  used  for  performing  this  calculation  is  listed  in  Appendix  A. 


3.3.  RESULTS  AND  DISCUSSION  FOR  ELECTRIC  FIELD  NORMAL  TO 
SURFACE 

Current  densities  of  escaping  electrons,  computed  as  described  in  Sec.  3.2, 
are  shown  in  Table  2  and  Fig.  3.5.  Each  value  of  i  -  I/I  was  calculated  using 
80x80x40  orbits,  whose  initial  velocity  components  u  ,u  ,  and  u  were 
chosen  as  described  in  Sec.  3.2,  and  with  points  on  the  orbits  calculated  at 
intervals  Ac  =  rr/45.  With  these  increments,  the  results  have  a  numerical 
accuracy  of  about  0.2%  or  better.  For  8  values  of  e  and  1 1  values  of  0,  the 
resulting  calculation  took  about  100  hr  total  on  a  Hewlett-Packard  1000F 
computer  with  Vector  Instruction  Set.  The  result  for  e  =  0  is  just  the 
analytic  result  i  =  cos  0.  To  see  why  this  is  so,  we  consider  the  electron 
orbit  shown  in  Fig.  3.6,  which  has  been  fictitiously  extended  so  as  to  pass 
through  the  surface  and  re-emerge  from  it.  In  the  absence  of  an  electric  field 
(e  =  0),  this  orbit  has  the  same  speed  at  the  re-emergence  point  C  as  at  the 
emission  point  A.  Since  we  have  also  assumed  that  the  emitted  velocity 
distribution  is  isotropic,  and  therefore  a  function  of  speed  only,  the  real  orbit, 
for  which  C  is  the  emission  point,  must  carry  the  same  population  as  would 
the  fictitious  re-emerged  orbit.  The  flux  crossing  the  reference  surface  DE, 
which  is  i  B,  is  therefore  the  same  as  if  such  passages  and  re-emergences 
actually  occurred,  and  is  the  same  as  if  another  reference  surface  FG,  also  1 
B,  were  emitting  electrons  having  the  same  velocity  distribution.  However,  in 
reality,  the  electrons  come  from  the  real  surface  HJ,  which  is  not  1  B,  and  all 
the  electron-orbit  guiding  centers  which  are  inside  any  given  magnetic-flux 
tube  through  OF.  will  also  be  inside  the  projection  of  the  same  flux  tube  onto 
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HJ,  and  the  ratio  of  the  intersection  areas  of  this  tube  with  HJ  and  DE  is  just 
sec  0.  The  ratio  of  escaping  to  emitted  flux  must  therefore  be  the  reciprocal 
of  this,  or  cos  0,  as  stated  above. 

Also  evident  in  Fig.  3.5  is  the  fact,  mentioned  in  Sec.  3.2,  that  when  e  is 
large  enough,  electron  escape  becomes  essentially  complete  except  when  0  is 
very  nearly  90°.  This  means  that  in  Shuttle  high-voltage  charging  conditions, 
for  which  30  ^  e  ^  160  (Sec.  3.2),  the  occurrence  of  high-voltage  charging  in 
marginal  circumstances  may  depend  very  strongly  on  the  precise  orientation  of 
a  surface.  A  slowly-rotating  surface  which  passes  through  tangentiality  to  B 
may  experience  a  sudden,  brief  high-voltage  charging  event.  For  the  same 
reason,  attempts  to  predict  high-voltage  charging  may  be  afflicted  by  "sensi¬ 
tivity  effects",  which  are  defined  as  large  changes  in  predicted  results  result¬ 
ing  from  small  changes  in  physical  input  parameters.  An  important  conse¬ 
quence  of  sensitivity  effects  here  is  that  if  one  attempts  to  predict  "worst- 
case"  charging  by  assuming  that  emitted  electrons  do  not  escape,  then  the 
resulting  predictions  are  likely  to  be  overly  pessimistic  most  of  the  time.  On 
the  other  hand,  if  one  assumes  that  they  do  escape,  correct  predictions  will  be 
obtained  almost  all  of  the  time,  but  occasionally  a  large  underestimate  of 
charging  will  occur.  A  different  kind  of  sensitivity  effect,  involving  ambient 
electron  distributions  which  are  almost  Maxwellians  having  temperatures 
close  to  the  "threshold  temperature  for  high-voltage  charging"  of  the  spacecraft 
surface  material,  was  identified  by  Laframboise  et  al  (1982a)  and  Lafram- 
boise  and  Kamitsuma  (1983);  see  also  Schnuelle  et  al  (1981,  Fig.  4). 
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The  results  for  i  in  Table  2  are  approximated  with  an  absolute  error  of 
2.5%  or  less  by  the  empirical  formula: 

f 

bj  =  i  +  i . 35c1  • 1394  exp  -  0.083725  fl  +tanh  [1.9732  in  q"~T~3  )]| 

] 

-0.07825  in  [1  +  (e/8 . 5) 1  * 78 1 4B]  i: 

j 

b2  =  0.38033e°-95892  exp|2.Q988|l  +  tanh  [1.49  in  (3^5)]}  ; 

c  =  in  (90°/ 0) ;  (3.11) 

i  =  cos  [90°exp(-b^c-b2C2)]. 

This  formula  also  has  the  correct  limiting  behavior  when  e  -►  0  or  00,  or  0  -* 
0°  or  90°.  In  order  to  calculate  the  escaping  flux  using  (3.11),  one  also  needs 
to  know  the  emitted  flux  I  .  For  Dhotoemission,  values  of  I  for  various  space- 
craft  materials  have  been  given  by  Feuerbacher  and  Fitton  (1972)  and  Grard 
(1973).  For  secondary  and  backscattered  emission,  analytic  approximations 
for  fractional  yields  (IQ/incident  flux)  as  functions  of  incident  energy  and 
direction  have  been  given  by  Laframboise  et  al  (1982a)  and  Laframboise  and 
Kamitsuma  (1983);  see  also  Katz  et  al  (1986).  The  results  presented  in  this 
Section  have  already  been  presented  by  Laframboise  (1985). 

In  a  real  situation,  F.  would  not  be  uniform,  but  would  decrease  with  dis¬ 
tance  from  the  surface,  contrary  to  our  assumptions.  Our  results  can  there¬ 
fore  be  expected  to  overestimate  electron  escape.  This  would  probably  net  be 
a  large  effect.  For  example,  the  spherical -sheath  thicknesses  mentioned  in 


Sec.  3.2  for  -lkV  and  -5kV  surface  potentials  are  2.6  and  6.1  m,  so 
variations  of  E  will  be  small  within  the  13  cm  gyroradius  of  a  3eV  emitted 
electron.  If  desired,  one  can  make  an  approximate  compensation  for  this 
nonuniformity  by  calculating  e  using  an  electric  field  value  which  is  averaged 
over  the  first  average-gyroradius  distance  from  the  surface. 


3.4.  THEORY  FOR  NON-NORMAL  ELECTRIC-FIELD  DIRECTIONS 

If  a  spacecraft  surface  is  charged  to  a  nonuniform  potential,  the  electric  field 
E  outside  it  has  a  nonzero  component  tangential  to  the  surface.  The  resulting 
problem  geometry  is  more  complicated,  and  two  additional  angles,  a  and  \py 
which  define  the  direction  of  the  electric  force  -eE,  need  to  be  specified 
(Fig.  3.7).  We  still  assume  that  E  and  Bare  uniform,  as  in  Sec.  3.2. 

It  may  then  happen  that  even  though  -eE  is  still  directed  away  from  the 
surface,  its  projection  along  the  direction  of  B  is  directed  toward  the  surface, 
as  shown  in  Fig.  3.8a.  This  happens  when  the  angle  between  -eE  and  either  B 
or  B,  whichever  is  outward,  is  greater  than  90°.  The  resulting  situation  is 
analogous  to  that  for  a  sailboat  tacking  against  the  wind  (Fig.  3.8b). 

When  this  situation  exists,  every  emitted  electron  returns  to  the  surface, 
and  the  escaping  flux  is  zero.  The  resulting  situation  then  differs  rrorn  that 
shown  in  Fig.  3.5  in  at  least  one  important  respect.  The  range  of  surface 
orientations  (relative  to  B)  for  which  electron  escape  is  entirely  prevented,  is 
no  longer  infinitesmal  (at  0  90  ),  but  finite,  and  this  may  greatly  enlarge 
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the  portion  of  a  spacecraft’s  surface  for  which  secondary-electron  or  photo¬ 
electron  escape  is  not  available  as  a  discharge  mechanism. 

Even  though  no  electrons  escape  under  these  conditions,  the  possibility 
exists  that  they  may  travel  a  relatively  long  distance,  equal  to  many  gyre- 
radii,  parallel  to  the  surface  before  returning.  This  may  produce  relatively 
large  surface  currents,  and  these  may  modify  substantially  the  charge  distrib¬ 
ution  on  the  spacecraft.  The  question  of  surface  currents  resulting  from  migr¬ 
ation  of  reimpacting  electrons  along  the  surface  remains  to  be  examined  in 
detail;  this  situation  is  discussed  in  S^  .-tion  6.  Here,  we  point  out  only  one 
general  feature  of  this  charge  migration.  The  general  motion  of  an  emitted 
electron  is  a  superposition  of  gyromotion  about  B,  an  E  x  B  drift,  and  an 
acceleration  in  the  B  or  -B  direction.  The  ExB  drift  direction  may  have 
either  an  inward  or  outward  normal  component.  If  it  is  outward,  migration 
distances  over  the  surface  will  be  larger,  and  3u  therefore  will  surface  cur¬ 
rents.  Since  B  is  in  the  (x,z)  plane,  the  z-component  of  E  x  B/B2  is 
-Bx  Ey/B2.  For  B^  >  0,  this  component  has  a  sign  opposite  to  that  of  E^. 
"Therefore,  when  E^  <  0  [region  (3)  in  Fig.  3.9],  the  migration  distances  of 
reimpacting  electrons,  and  therefore  also  the  surface  currents  produced  by 
them,  are  likely  to  be  much  larger  than  when  E^  >  0  [region  (2)  in  Fig.  3.9]. 
These  surface  currents  will  be  primarily  along  the  tangential  projection  of  the 
FI  x  B  drift  d  ireeticn,  rather  than  of  -eF.,  so  they  may  have  little  effect  on 
helping  to  discharge  the  spacecraft.  For  some  orientations  of  E  and  B,  emitted 
electrons  may  gain  enough  kinetic  energy  before  reimpacting  to  cause  substan¬ 
tial  amounts  of  secondary  emission  to  occur  when  they  do  so.  Surface  currents 
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then  may  be  greatly  increased  by  electron  multiplier"  effects  (Goodrich  and 
Wiley,  1961;  Wiley  and  Hendee,  1962). 


Our  predictions  are  based  on  the  assumption  that  F.  and  B  are  spatially 
uniform  and  time-independent.  However,  if  average  migration  distances 
become  equa.  to  many  gyroradii,  it  is  then  more  likely  that  spatial  uniformity 
will  be  seriously  violated,  and  many  of  our  predictions,  especially  the  quanti¬ 
tative  ones  presented  m  Sec.  3.5,  may  then  become  unreliable.  The  same 
circumstances  increase  the  possibility  that  time-dependent  phenomena  (plasma 
turbulence)  may  also  affect  our  predictions.  This  possibility  may  be  remote, 
because  the  spontaneous  plasma  oscillations  observed  recently  in  the  disturbed 
region  around  the  Shuttle  are  of  relatively  small  amplitude.  They  involve 
relative  density  fluctuations  (An/n)  of  at  most  a  few  percent,  and  the  largest 
amplitudes  of  these  fluctuations  generally  do  not  occur  close  to  the  space¬ 
craft’s  surfaces  (Murphy,  1935;  Murphy  et  al,  1986).  However,  this  question 
remains  to  be  thoroughly  explored. 

Before  we  present  computed  results,  we  need  to  determine,  in  terms  of  the 
angles  6,  a,  and  ^  in  Fig.  3.7,  when  the  projection,  of  -eF.  along  B  is  directed 
toward  the  surface.  This  projection  is: 

(eE»B)B/B2  (3.12) 

The  z-comnonent  of  this  is: 

i 

-e(E  B  +  E  B  )  B  /B2 

x  y  z  z  z 

=  -  eE  co mO  cos  or  (1  +  tan  a  cos  \p  tan  0)  (3. 13) 
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(3.14) 


Since  -eE  cos2  0  cos  a  >  0,  this  means  that  escape  is  prevented  if: 
tan  a  cos  p  tan  0  <  1 , 

or  eauivalentiv: 

I  J 

0  <  tan'1  (-cot  a  sec  p) ,  for  0  T  p  <(  90°; 

0  9  tan  1  (-cot  a  sec  p) ,  for  90°  <  y  ^  180°. 


(3.15) 


3.5.  RESULTS  AND  DISCUSSION  FOR  NON-NORMAL  F.LF.CTRIC-FIELD 
DIRECTIONS 

Figures  3.10  -  3.17  show  escaping  electron  current  densities  i  =  i (a,p,0,e) 
when  -eE  is  not  normal  to  the  spacecraft  surface  (or  ^  0°).  Details  of  the 
computations  of  these  results  are  the  same  as  those  given  in  Sections  3.2  and 
3.3,  except  that  64  x  64  x  32  orbits  were  used  for  calculating  each  value 
of  i,  and  the  results  have  a  numerical  accuracy  of  about  0.4%  or  better. 

In  Figs.  3. 10  and  3.11,  p  “  0°,  so  -eE,  B,  and  the  surface  normal  are 
coplanar.  The  E  x  B  drift  direction  is  therefore  tangential  to  the  surface.  In 

Fig.  3.10,  no  electrons  escape  (i  -  0)  when  -90°  ^  0  <  -60°  because  the 
acceleration  of  all  electrons  along  B  is  toward  the  surface,  as  discussed  in 
Sec.  3.4.  For  0  >  -60°,  this  acceleration  is  away  from  the  surface,  but  it  is 
larger  for  0  >  0°  than  for  0  <  0°,  so  increasing  the  electric  field  magnitude  e 
increases  electron  escape  more  for  0  >  0°.  For  0  >  -  60°,  l  -  cos  0  when  c  - 
0,  as  was  the  case  for  -eE  normal  to  the  surface  (Fig.  3.5).  Some  care  is 
needed  in  defining  what  is  meant  by  the  case  "c  ~  0".  For  90°  ^  0  <  -60°, 
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i  =  0  in  the  linn, it  e  ►  0  f  ,  but  in  the  limit  e  ►  ()-,the  electric  field  is  reversed, 
so  i  =  cos  0  in  tins  limit,  and  i  therefore  has  a  discontinuity  at  e  -  0.  The 
discontinuity  is  reversed  for  0  >  60°.  However,  we  are  interested  here 
primarily  in  cases  when  e  /  U  (negat  iveiy- charged  spacecraft  surfaces),  so  in 
this  work,  we  take  T  -  0"  to  mean  the  limiting  case  t:  -►  0+.  In  Fig.  3.11,  a 
has  been  increased  from  30°  to  60°,  so  the  effects  just  discussed  in  connec¬ 
tion  with  Pig.  3.10  are  seen  again,  but  more  strongly.  This  time,  escape  is 
suppressed  completely  for  -90°  ^  0  <  -30°. 


In  Fig.  3.12,  the  acceleration  of  electrons  along  B  is  toward  the  surface, 
and  therefore  i  =  0,  for  0  <  tan'1  (-VS)  =  -67.79°,  as  given  by  Eq.  (3.15). 
Also,  we  now  have  gf  0°,  so  the  electric  force  vector  is  no  longer  in  the 
same  plane  as  B  and  the  surface  normal.  As  a  result,  the  P  xR  drift  now  has 
a  nonzero  normal  component.  This  decreases  electron  escape  for  0  <  0°,  and 
increases  it  for  0  >  0°.  It  also  causes  the  escape  to  remain  nonzero  at  0  = 
90°.  As  before,  the  larger  outward  acceleration  along  B  aiso  increases 
escape  for  0  >  0  ’.  For  0  just  larger  than  67.79°,  we  see  that  escape  is 
suppre53sed  almost  complete’y  for  larger  values  of  c;  this  is  because  the 
inward  direction  of  the  E  x  B  normal  comoonent  causes  most  electrons  to 


reimpact  the  surface  during  the  first,  gyroperiod  after  emission.  In  Fig. 

3.13,  a  has  been  increased  from  30°  to  60°,  with  consequent  enhancement  of 


the  effects  just  discussed.  We  now  have  i*-0  for  0  <  tdn“'( 


OQ  OO0 


T.n  Figs.  3.  Id  and  3.15,  y  --  30°,  and  the  projection  cf  eE  along  B  is  away 
from  the  surface  for  ail  0,  so  suppression  of  electron  escape  by  deceleration 
along  B  does  not  occur.  However,  for  0  close  to  30°,  the  effect  of  the  inward 
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direction  of  the  E  x  B  norma:  component  overcomes  the  effect  of  the  outward 
acceleration  along  B,  especially  because  B  is  now  nearly  tangential  to  the 
surface,  and  therefore  electron  escape  is  effectively  suppressed  for  larger 
values  of  e.  In  Fig.  3.15,  suppression  for  c  -  20  is  essentially  complete  over 
a  range  of  0  values  extending  mere  than  30°  on  either  side  of  0  =  -90°. 

In  Figs.  3.16  and  3.  17,  we  have  ip  =  135°,  and  the  effects  of  electron 
deceleration  along  B  and  of  El  x  B  drift  now  suppress  electron  escape  at 
opposite  ends  of  the  range  of  0.  From  Fa.  (3.15),  we  now  have  i  =  0  for  0  )> 
tan_1(V6)  =  67.79°  and  G  >  tan' 1  ( V2 / 3 )  =  39.23°,  respectively.  In  Fig.  3.16, 
the  E  x  B  effect  significantly  enhances  electron  escape  as  0  increases,  just 
before  the  deceleration  effect  cuts  it  ofF 

The  results  shown  in  Fig.  3.5  and  Figs.  3.10-3.17  encompass,  albeit 
rather  sparsely,  the  entire  range  of  possible  directions  of  B  and  E  for  a  up  to 
60°.  To  see  this,  we  first  note  that  in  the  important  case  where  -eE  is  normal 
to  the  surface,  the  i  values  for  -90°  ^  0  s  0°  can  be  generated  from  those  for 
tne  range  0  ^  G  ^  90  ,  which  is  covered  in  Fig.  3.5;  this  can  be  seen  by 
rotating  the  B  vector  in  Fig.  3.7  by  180°  about  the  z  axis.  Secondly,  the  i 
values  for  a  =  30°  and  60°  and  =  180°,  225°,  270°,  and  315°  can  be 
generated  from  those  in  Figs.  3.10-3.17  by  rotating  the  B  and  -eE  vectors  ir. 
Fig.  3.7  together  about  the  z  axis  by  180°.  The  effect  of  this  is  to  increase 
all  the  \p  values  by  180°  and  also  to  reverse  the  sign  of  0.  Finally,  we  can 
obtain  the  :  values  for  cases  where  B  is  reversed  by  first  noting  that  reversal 
of  B  implies  reversal  of  both  the  £  coordinate  in  Fig.  3.7  (in  order-  that  £ 
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remain  parallel  to  B)  and  the  i)  coordinate  [in  order'  that  the  (£,y,q)  axes 


remain  right-handed].  The  quantities  sin  0,  ccs  0,  uq^.,  uq^,  e^?  and  in 
£q.  (3.8)  will  then  all  reverse.  To  keep  z  =  z(r)  in  Fq.  (3.8)  unchanged,  we 

require  also  that  n  reverse,  and  we  therefore  read  re  that  e  and  u  also  be 
t  ■  y  oy 

reversed.  However,  the  emitted  velocity  distribution  is  symmetric  in  u  , 
and  reversing  involves  replacing  the  angle  \p  by  360°  -  \p,  and  this 
replacement  gives  hack  the  same  set  of  \p  values  for  which  our  computations 

already  give  i. 


We  can  summarize  the  results  in  Figs.  3. 10-3.17  bv  noting  that  when  a  t 
0.  two  new  mechanisms,  which  were  not  present  when  -eE  was  normal  to  the 
surface  (Fig.  3.5),  can  suppress  electron  escape.  These  are:  an  inward 
normal  component  of  E  x  B,  if  e  is  strong  enough,  and  a  decelerating 
projection  of  eE  along  B,  for  any  e  >  0.  These  may  act  at  the  same  end  or  at 
opposite  ends  of  the  range  of  magnetic-field  directions  -90  ^  0  £  90  .  These 
mechanisms  can  greatly  enlarge  the  range  of  surface  orientations  for  which 
escape  is  suppressed. 
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3.6.  CM  .GUI  ATION  OF  ESCAP!N(VEI  .KCTRON  DENSITIES 


Once  the  fluxes  of  escaping  electrons  are  known  (Sections  3. 3-3. 5),  a 
simple,  inexpensive,  approximate  calculation  of  their  space-charge  density 
distribution  can  be  set  up.  The  proposed  method  is  as  follows:  (1)  ignore  the 
gyromotion  of  electrons  once  they  have  escaped.  Their  motion  then  involves: 

(a)  an  acceleration  along  magnetic  field  lines,  of  amount  (e/m)E»B/B  (b)  a 
drift  motion  of  velocity  E  x  B/B2  across  magnetic  field  lines.  (2)  Integrate 
enough  of  the  trajectories  defined  by  this  motion  (i.e.  their  guiding-center 
trajectories)  to  define  trajectory  tubes  whose  cross-section  at  any  point  can  be 
calculated  with  sufficient  accuracy;  the  method  described  by  Laframboise  et 
ai  (1982b,  Sec.  /),  can  be  used  to  calculate  the  area  of  a  trajectory  tube 
without  reference  to  neighbouring  trajectories.  (3)  Calculate  their  space- 
charge  density  n(r)  at  any  point  by  (a)  ignoring  the  "thermal"  spread  of  their 
velocities  (b)  then  invoking  the  fact  that  their  density  x  their  velocity  [as 
given  by  the  orbit  integration  mentioned  in  (2)],  x  the  cross-sectional  area 
A(r  )  of  the  trajectory  tube  (which  must  be  calculated  in  a  plane  1  the  trajec¬ 
tory)  at  the  point  r  in  question,  =  a  constant  (whose  value  is  given  by  the 
initial  conditions  at  the  point  on  the  spacecraft  where  the  trajectory  origi¬ 
nates)  (c)  finding  their  velocity  at  the  point  in  question  by  using  energy  con¬ 
servation,  together  with  the  values  of  electric  potential  c p(r)  and  <pQ  at  that 
point  and  the  emission  point,  and  their  assumed  velocity  vq  at  the  emission 
point.  The  result  is: 


n(r)  =  n  v  A  / 
o  o  o 


A(r ) 


2  +  (2e/rn)  [$(r  ) 


(3.16) 
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where  nQv0  is  the  escaping  flux  calculated  in  Sections  3. 3-3. 5.  At  most 
positions,  n(r)  will  be  insensitive  to  the  precise  value  assumed  for  vq2; 

assuming  that  vq  =  the  one-sided  thermal  speed  (SkT/irm'^  will  suffice  for 
most  purposes. 


37 


u . 


HIGH-VOLT  AGP.  CHARGING  OF  A  DIELECTRIC -DOVER  ED 
SPACECRAFT:  A  WAKE- INDUCED  BARRIER -EFFECT  MECHANISM. 


4.1. 


INTRODUCTION  TO  THE  WAKE  INF'  ICED- BARRIER-EFFECT 
PROBLEM 


In  Section  3,  we  examined  the  suppression  of  emitted-electron  escape  by 
magnetic  fields.  On  a  dielectric-covered  spacecraft,  this  suppression  can  be 
expected  to  result  ir.  local  high-voltage  charging  in  certain  kinds  of  auroral- 
plasma  conditions  which  would  not  produce  overall  charging.  This  point  is 
discussed  further  in  Sections  1  and  6. 

In  this  Section,  we  examine  a  different  effect,  namely  the  wake-induced- 
barrier  effect,  which  produces  a  similar  result,  i.e.  suppression  of  the  escape 
of  emitted  electrons.  Our  work  here  is  directed  specifically  to  the  situation  of  a 
large  spacecraft  with  dielectric  exterior  surfaces  either  in  the  auroral  plasma  or 
emitting  an  electron  beam  as  in  the  SEP  AC  Spaeelab-1  experiment. 

The  study  of  Parks  and  Katz  (1981),  Katz  and  Parks  (19831  concerned  overall 
high-voltage  charging  of  a  conductive  spacecraft.  Our  work  in  this  Section 
snows  that  local  high-voltage  charging  on  a  large  dielectric  spacecraft  can  be 
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expected  to  occur  in  the  auroral  ionosphere  in  conditions  different  than  those 
already  known  to  have  produced  high-voltage  charging  on  other  spacecraft.  In 
Section  4.2,  we  describe  the  wake-induced-barrier  effect,  including  the 
differences  between  it  and  the  barrier  effect  on  geosynchronous-altitude 
satellites.  In  Section  4.3,  we  develop  theory  for  calculating  spacecraft 
surface  potentials  produced  by  the  wake-induced-barrier  effect.  In  Sections  4.4 
and  4.5,  we  apply  this  theory  to  spacecraft  in  the  auroral  plasma  or  emitting  an 
electron  beam,  respectively. 
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4.2.  DESCRIPTION  OF  THE  WAKE-INDUCED-BARRIER-EFFECT 
MECHANISM 

Ir.  low-orbit  conditions,  the  plasma  environment  of  a  spacecraft  is  usually 
''mesothermal",  i.e.  it  has  ion  mean  thermal  speed  <<  spacecraft  speed  <X 
electron  mean  thermal  speed.  Under  these  conditions,  the  ambient  ion  flux  in  the 
soacecraft  reference  frame  is  mostly  from  the  upstream  direction,  whereas 
ambient  electrons  come  morc-or-less  equally  from  all  directions  [unless  the 
ambient  electron  distribution  itself  is  highly  anisotropic;  see  Fennell  et  al 
(1981).  Lin  and  Hoffman  (1982),  and  Yeh  and  Gussenhoven  (1987)].  Relatively 
few  ions  are  then  able  to  reach  the  region  just  downstream  of  the  spacecraft,  and 
a  “wake"  region  of  unbalanced  electron  space  charge  forms  there.  This  wake 
region  is  more-or-iess  "wedge-shaped"  (or  cone-shaped,  depending  on  the 
spacecraft  geometry).  If  the  ratio  of  spacecraft  size  to  (electron)  Debye  length 
is  large  enough,  this  unbalanced  space  charge  is  believed  to  produce  a  potential 
minimum  downstream  of  the  spacecraft.  The  existence  of  such  a  minimum  has 
been  predicted  by  many  authors  including  Maslennikov  and  Sigov  (1965,  1967, 
1969),  Maslennikov  et  al  (1968),  Taylor  (1967),  Kiel  et  al  (1968),  Call  (1969, 
Figs.  3. 1-3.3.  3.14,  3.21.  and  3.23),  Gurevich  et  al  (1970),  and  Vaglio- 
Laurin  and  Miller  (1970).  These  papers  are  part  of  an  extensive  literature  on 
wakes  of  spacecraft  in  low  Earth  orbit.  Reviews  of  this  literature  have  been 


40 


given  by  Brundir,  (1963),  Gurevich  et  al  (1970),  Samir  and  Stone  (1986),  and 
Stone  and  Samir  (1986). 

In  situations  studied  by  the  above-mentioned  authors,  the  spacecraft  is 
assumed  to  be  conductive,  and  therefore  its  potential  relative  to  space  will 
"float"  at  a  uniform  value  determined  by  total,  rather  than  local,  current  balance. 
If  the  ambient  electrons  are  cold  (kTg  -  0.1  eV),  then  this  floating  potential 
will  be  a  few  tenths  of  a  volt  negative  with  respect  to  space.  The  resulting 
potential  distribution  downstream  of  the  spacecraft  has  the  general  appearance 
shown  in  Fig.  4.1(a). 

We  now  consider  situations  which  will  arise  during  polar-orbit  flights  planned 
for  the  Shuttle  Orbiter.  Such  situations  will  differ  from  the  one  just  described 
in  two  important  ways.  First,  the  Orbiter’s  surface  is  mostly  dielectric  rather 
than  conductive.  Second,  the  auroral -rone  ambient  plasma  often  contains  a  high- 
energy  "hot"  electron  constituent  whose  equivalent  kinetic  temperature  is 
generally  in  the  range  1-15  keV  (Gussenhoven  et  al,  1985;  Yeh  and 
G  ussenho  ven ,  1987). 

In  such  situations,  each  small  dielectric  portion  of  the  spacecraft’s  total 
surface  will  float  separately  at  a  potential  determined  by  local  rather  tnan  total 


current  balance.  If  enough  hot  auroral  electrons  are  present,  the  situation  shown 
in  Fig.  4.1(a)  will  no  longer  he  able  to  maintain  itself.  To  see  why,  we  note  that 
the  negative  potentials  present  in  the  wake  will  exclude  the  cold  electron 
population  from  it  almost  completely,  so  the  hot  electrons  will  form  almost  its 
entire  population.  The  potential  at  the  minimum  shown  in  Fig.  4.1(a)  will  then 
be  more  negative  than  the  potentials  on  adjacent  surfaces  by  an  amount  which  is 
much  larger  than  the  average  emission  energies  of  either  secondary  electrons  or 
photoelectrons.  The  resulting  potential  barrier  will  therefore  reflect  these 
electrons  and  cause  them  to  return  to  these  surfaces  at  locations  generally  close 
to  then  emission  points.  As  a  result,  the  secondary-emission  or  photoemission 
current  contributions  will  become  unavailable  as  a  means  of  discharging  these 
surfaces,  which  will  then  charge  progressively  more  negatively  until  surface 
potentials  "just  inside"  the  potential  minimum  have  become  almost  as  negative  as 
at  the  minimum.  A  finite  fraction  of  the  emitted  secondaries  and  photoelectrons 
will  now  escape,  and  these  surfaces  will  now  arrive  at  a  steady-state  "floating" 
condition.  We  shall  see  later  (Sec.  4.4)  that  on  a  large  enough  spacecraft, 
surface  potentials  more  than  100  volts  negative  with  respect  to  space  can  be 
expected  because  of  this  mechanism.  In  the  remainder  of  this  Section,  we 
discuss  some  other  features  of  it. 
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One  such  feature,  which  is  characteristic  of  it,  is  that  the  location  of  the 
potential  minimum  will  migrate  inward  until  it  almost  reaches  the  spacecraft. 
The  general  appearance  of  the  equipcter.tial  surfaces  downstream  of  the 
spacecraft  will  then  be  as  shown  in  Fig.  4.1b. 

A  similar  phenomenon,  which  often  controls  the  high-voltage  differential 
charging  of  geosynchronous-altitude  satellites,  is  called  the  "barrie*-  effect”. 

This  effect  was  first  predicted  by  Fa'nleson  (1973)  and  observed  by  Whipple 
(1976).  It  was  then  rediscovered  independently  in  a  numerical  simulation  by 
Katz  et  al  (1979).  It  has  also  been  studied  by  numerous  other  authors 
(Prokopenko  and  Laframboise,  1977,  1980;  Laframboise  and  Prokopenko,  '977; 
Besse  and  Rubin,  1980;  Olsen  et  al,  1981;  Purvis,  1982;  Katz  and  .Vandell, 
1982;  Laframboise  et  al,  1982;  Olsen  and  Purvis,  1983).  However,  two 
features  of  the  present  situation  are  very  different  than  in  the  geosynchronous- 
altitude  situation.  In  that  situation,  the  Debye  length  of  the  ambient  plasma  was 
much  larger  than  the  spacecraft,  and  therefore  the  barrier  could  not  be  produced 
by  space-charge  effects.  Barrier  formation  then  occurred,  over  some  part  of  the 
spacecraft  surface,  because  another  part  of  it  had  become  charged  to  a  much 
larger  negative  potential.  For  this  to  occur,  the  other  part  had  to  have  a 
secondary-electron  yield  less  than  unity  (for  the  existing  incident  electron 
distribution)  and  also  had  to  be  shaded,  so  that  neither  secondary-electron 
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emission  nor  photoemission  could  discharge  it.  In  our  situation,  there  is  no 
requirement  that  a  larger-potential  surface  exist  somewhere  else,  and  the 
barrier  effect,  instead  of  limiting  the  differential  charging,  actually  produces  it. 
Secondly,  in  the  geosynchronous-altitude  situation,  the  barrier  effect  is 
characterized  by  a  saddle  point  in  the  potential  outside  the  spacecraft,  and  as  the 
surfaces  inside  the  barrier  charge  more  and  more  negatively,  the  saddle  point 
migrates  closer  and  closer  to  the  spacecraft  until  it  almost  reaches  it,  and  the 
potential  difference  between  the  saddle  point  and  the  surface  adjacent  to  it 
decreases  to  a  few  volts,  permitting  a  significant  fraction  of  the  emitted 
electrons  to  escape.  In  our  situation,  there  is  no  saddle  point.  Instead,  as 
mentioned  above,  there  is  a  minimum  in  the  potential  outside  the  spacecraft,  and 
as  the  surfaces  inside  the  barrier  charge  more  and  more  negatively,  this 
minimum  migrates  almost  to  the  surface  [Fig.  4.1(b)].  More  specifically,  the 
entire  surface-potential  distribution  must  adjust  itself  until  a  "ridge",  having  an 
almost  uniform  (negative)  "height"  comparable  to  the  emitted-eleetron  mean 
energy,  has  formed  outside  the  surface.  Since  this  mean  energy  is  only  a  few 
eV,  this  "ridge"  will  form  close  to  the  surface  [Fig.  4.1(b)].  We  see  that  a  good 
estimate  of  the  largest  (negative)  value  of  surface  potential  is  given  simply  by 
the  largest  value  of  potential  in  the  wake. 
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We  now  develop  theory  for'  calculating  what  values  of  surface  potential  will  be 
achieved  by  this  mechanism,  as  a  function  of  spacecraft  size  and  ambient  plasma 
conditions.  This  is  done  in  the  next  Section. 
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4.3.  SURFACE-POTENTIAL  CALCULATION 


We  wish  to  develop  a  method  of  estimating  the  largest  (negative)  value  of 
surface  potential  relative  to  space,  as  a  function  of  position  on  the  spacecraft’s 
wake-region  surfaces.  In  order  to  do  this,  we  make  a  set  of  approximations 
which  permit  relatively  simple,  partly  analytic  solutions  for  this  potential,  while 
retaining  the  most  important  physical  features  of  the  problem.  Such  a 
procedure  is  especially  valuable  in  view  of  the  complexity  of  an  exact  treatment 
(Section  4.4;  Section  5).  These  approximations  are  as  follows. 

As  mentioned  in  Section  4.2,  we  assume  that  the  most  negative  potential  in 
the  wake  region  just  downstream  of  the  spacecraft  is  a  good  approximation  to  the 
most  negative  potential  on  the  spacecraft. 

We  assume  that  the  spacecraft  is  an  infinite  cylinder.  However,  for  most  of 

our  purposes,  its  detailed  shape  will  not  be  important;  its  most  important 

dimension  will  be  its  least  dimension  transverse  to  its  motion.  We  assume  that 

its  surface  is  made  of  dielectric  material,  and  it  is  immersed  in  a  collisionless 

plasma  crossflow  which  has  a  transverse  drift  velocity  U  =  8  km/sec  relative  to 

the  spacecraft,  corresponding  to  the  speed  of  low  Earth  orbit.  The  ions  in  the 

plasma  are  assumed  to  be  0+  at  a  temperature  of  0.1  eV.  Their  most-probable 

thermal  speed  (2kT./mJ^  is  then  1.09  km/sec,  and  this  is  much  less  than  U, 

l  i 
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so  we  can  expect  the  wake  to  extend  for  many  spacecraft  diameters  in  the 
downstream  direction  (Fig.  4.1a).  The  electrons  in  the  plasma  are  assumed  to 
include  a  "coid,,  (0.1  eV)  and  a  “hot"  (1  keV  or  more)  Maxwellian  constituent, 
representing  auroral  plasma  conditions.  This  is  a  very  rough  approximation  to 
observed  auroral -elec iron  velocity  distributions,  which  can  be  both  very  non- 
Maxwell  lan  and  very  anisotropic  (Sec.  4.2),  but  it  is  adequate  for  our  purposes. 

We  ignore  magnetic-field  effects  on  charged-particle  motions.  This 
approximation  should  be  acceptable  for  our  purposes  because  the  average 
gyroradii  of  the  hot  electrons  and  the  ions  are  a  few  meters  or  larger,  especially 
in  a  high-voltage  wake  region,  and  the  cold  electrons  are  almost  completely 
excluded  from  the  same  region.  Furthermore,  collection  of  cold  electrons  by  a 
negatively-charged  spacecraft  is  very  small,  and  therefore  their  density  is  well- 
approximated  by  a  Boltzmann  factor,  independently  of  magnetic-field  effects. 

We  assume  that  the  r  otation  rate  of  the  spacecraft  is  slow  enough  that 
transfer  of  surface  charge  into  the  spacecraft’s  wake  region  by  rotation  ("spin 
charging";  Parks  and  Katz,  1983)  is  negligible. 

Our  most  severe  approximations  concern  the  ion-density  distribution  in  the 
spacecraft  wake.  Exact  calculations  of  these  distributions  are  very  difficult 
(Sec.  4.4;  Sec.  5)  We  assume  that  the  large-voltage  region  of  the  spacecraft’s 


47 


wake  is  sufficiently  elongated  that  variations  in  potential  parallel  to  the  wake 
axis  (the  x  axis  in  Figure  4.1)  are  negligible  in  comparison  with  variations 
perpendicular  to  the  same  axis.  This  allows  us  to  write  a  one-dimensional 
Poisson  equation  for  potentials  in  this  region.  We  further  assume  that  outside 
this  region,  the  potential  <p  is  equal  to  space  potential,  which  we  take  as  our 
zero,  and  the  ion  and  electron  number  densities  n.  and  n^  are  equal.  We  assume 
that  the  cold-electron  density  decreases  abruptly  to  zero  at  the  edges  of  this 
region,  so  that  the  only  electrons  inside  it  are  the  hot  ones.  We  ignore,  for 
now,  the  effect  of  electron  collection  by  the  spacecraft  on  the  space-charge 
density  of  the  hot  electrons,  and  we  also  ignore  the  density  contribution  made  by 
the  emitted  electrons,  whether  these  return  to  the  surface  or  escape.  We  make 
an  approximate  correction  for  these  two  effects  later.  With  our  other 
approximations,  these  imply  that  the  electron  density  in  this  region  is  given  by  a 
Boltzmann  factor.  We  assume  that  the  ion  density  in  the  same  region  is  either 
zero  (in  the  crudest  approximation)  or  has  a  uniform  value  which  we  estimate 
later  in  this  Section  using  a  model-potential  calculation  (in  the  next 
approximation) . 

With  these  approximations,  Poisson’s  equation  in  the  large-voltage  wake 
region  (the  "wake  potential  well")  becomes: 
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d20  _  n  ooehe 


d££  _ 
dy2 


exp 


e<p 


(4.1) 


where  y  is  distance  perpendicular  to  the  wake  centerline  (Fig.  4.1),  Te^  and 
r‘coeh  are  ternPerature  and  ambient  number  density  of  the  hot  electron 
constituent,  e  is  the  magnitude  of  unit  electronic  charge,  e0  is  the  permittivity 
of  space,  and  k  is  Boltzmann’s  constant.  Boundary  conditions  on  the  potential 
are: 


0  =  0  when  y  =  ±  y  (4.2) 

^  1  ;max 

where  ymax  is  the  half-width  of  the  wake  potential  well  just  downstream  of  the 
spacecraft.  Using  symmetry,  we  have: 

4>  =  4>0  <  0,  ^  =  0  when  y  =  0.  (4.3) 

By  our  assumptions,  y__v  is  just  the  spacecraft  radius  r  .  We  discuss  later 

m  9X  o 

how  to  obtain  an  "improved"  value  for  y 

r  7max 


We  define  the  following  dimensionless  variables: 

?  =  y/XDeh  =  y/te°kTeh/e‘naoeh>i 

(4.4! 

+  =  -ei£/kTeh 


where  Xq  ^  is  the  hot-electron  Debye  length.  Equations  (4.1)  -  (4.3)  then 
become: 
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d2  vl'/d£2  =  -  exp  (-44 

(4.5) 

*  =  0  when  ?  =  ±  fmax 

(4.6) 

^  =  'Fo  2  0,  d'F/df  =  0  when  £  =  0 

(4.7) 

Equation  (4.5)  is  nonlinear,  but  it  can  be  solved  by  multiplying  both  sides  by 
2  d^/ d£  and  noting  that  its  left  side  then  equals  (d/d£)  (d^/d£)2.  Its  solution 
subject  to  Eq.  (4.7)  then  is: 

£  =  ±  [2  expl^o)]^  Arctan  [expl^o  -  40  -1]^  (4.8) 

Using  (4.6)  now  gives  us,  in  implicit  form,  the  desired  relation  between  the 

normalized  wake  half-width  £  and  the  normalized  largest  negative  wake 

max 

potential  'I'q,  as  follows: 

6max  =  C2  expj^o)]^  Arctan  (exp^o)  -  1]^  (4.9) 

For  various  values  of  £  ,  the  solid  curves  labeled  "0"  in  Fig.  4.2  show  'I'  as  a 

function  of  the  normalized  transverse  distance  £  from  the  wake  centerline  in  the 
wake  potential  well,  as  given  by  Eq.  (4.8).  For  narrow  wakes  (small  £max  and 
4/0),  the  potential  variation  across  the  wake  is  nearly  parabolic,  but  for  wakes 
whose  width  is  many  hot-electron  Debye  lengths,  the  decrease  in  electron  density 
in  the  wake  produces  an  increase  in  local  Debye  length,  with  a  consequent 
"flattening"  of  the  potential  profile  in  the  well,  and  a  slowing-down  of  the  rate  of 
increase  of  well  depth  'I'0,  from  a  limiting  dependence  of  the  form: 
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*•  -  i  Tax2 


for  small  £ 


(4.10) 


-  ln<2  Tax*2*2’  for  larSe  Tax 


(4.11) 


The  dependence  of  'i'o  on  £max  given  by  Eq.  (4.9)  is  plotted  as  the  uppermost 
solid  curve,  labeled  "0",  in  Fig.  4.3.  Equations  (4.10)  and  (4.1  1)  are  plotted  in 
Fig.  4.3  as  dotted  curves.  Again,  'I'o  will  be  almost  equal  to  the  largest 
normalized  negative  surface  potential  0  on  the  spacecraft. 


Equation  (4.10)  is  equivalent  to: 

0o  ~ ►  -v  2  e  n  ,  /2e0 
y  max  ooeh  u 


(4.12) 


and  Eq.  (4.12)  indicates  that  for  potentials  small  compared  with  the  hot-electron 
temperature,  the  potential  produced  by  the  wake-induced-barrier  mechanism 
depends  only  on  the  ambient  density,  and  not  the  temperature,  of  the  hot 
electrons. 


The  most  important  approximation  in  our  calculation  so  far  is  the  omission  of 
ion  density  in  Eq.  (4.1).  We  now  estimate  the  effects  of  including  ion  density. 

If  nwl  is  ion  density  in  the  wake,  including  it  would  result  in  replacement  of  Eq. 
(4.5)  by: 
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(4.13) 


d2vT  _  ‘ '  wi 
d£2  _  n 


coeh 


If  n^.  were  uniform ,  it  could  never  be  large  enough  (for  a  given  'I')  to  reverse 


the  sign  of  d2vT/d£2,  because  the  solutions  of  (4.13)  would  then  become 
unphvsicai.  Therefore,  such  a  density  contribution  could  produce  only  a  further 
"flattening"  of  the  wake  potential  profiles,  and  a  limiting  of  ^  to  values  less  than 
that  at  which  this  sign-reversal  would  occur,  i.e: 


T'  <  ^  =  In  (n  , /n  .) 

n  ooeh  wi 


(4.14) 


We  have  solved  Eq.  (4.13)  numerically  for  various  values  of  ^max  and  n-wj/ncoef1- 
The  computer  program  used  for  doing  this  is  listed  in  Appendix  B.  Resulting 
wake  potential  profiles  are  shown  as  dotted  curves  in  Fig.  4.2.  Figure  4.3 
shows,  for  various  values  of  n^./n^^,  the  resulting  dependences  of  'I'o  on  £max- 
The  effects  noted  following  Eq.  (4.13)  are  clearly  evident  in  these  Figures. 
Equations  (4.11)  and  (4.14)  are  similar  to  Eqs.  (2)  and  (4a),  respectively,  of 
Vaglio-Laurin  and  Miller  (1970),  but  our  application  here  is  to  much  larger 
wake  potentials  than  theirs. 


We  now  introduce  an  approximate  correction  for  effects  on  electron  density 
arising  from  emitted-electron  space  charge  and  ambient-electron  collection  by 
the  spacecraft.  In  order  to  do  this,  we  ignore  magnetic-field  effects  on  electron 
motion  across  the  barrier,  and  we  assume  that  in  some  region  surrounding  a 
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point  on  the  barrier,  the  potential  distribution  has  planar  symmetry.  We  follow 
the  usual  procedure  of  approximating  the  emitted-electron  velocity  distribution, 
in  both  speed  and  direction,  by  a  Maxwellian  (Sec.  4.1).  The  velocity 
distributions  of  ambient  and  emitted  electrons  at  the  barrier  are  then  half- 
Maxwellians  directed  toward  and  away  from  the  spacecraft,  respectively.  If  we 
also  ignore  "hopping-conduction"  currents  of  emitted  electrons  along  the 
spacecraft  surfaces  inside  the  barrier,  then  the  inward  and  outward  fluxes 
carried  by  these  distributions  must  be  equal.  Since  the  emitted  electrons  have  a 
much  smaller  average  velocity  than  the  ambient  ones,  their  space-charge 
density  in  the  region  of  the  barrier  must  be  much  larger.  This  constitutes  a 
major  increase  in  wake-region  electron  density,  and  we  shall  see  later  (Sections 
4.4  and  4.5)  that  it  greatly  increases  predicted  wake-region  negative  surface 
potentials.  If  <p^  is  local  barrier  potential  and  Ts  is  the  temperature  of 
secondary  electrons  or  photoelectrons  emitted  from  the  surface,  then  with  the 
above-mentioned  approximations,  the  total  electron  density  at  locations  outside 
the  barrier  is: 


n  ((b)  =  n  , 
e  ooeh 


1 

7 


Teh'i 


exp 


e$b  e(cp-<pb) 

kT  + 


+  eXD  I 


exp  1  k’i  eh 


erfc 


‘  eh 


e(0-<^) 


erfc 


~v r 


kT 


eh 


(4.15) 
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where  erfc(x)  =  1  -  erf(x)  is  the  complementary  error  function.  At  the 
barrier,  Eq.  (4.15)  reduces  to: 

ne(V 


' 

-3° 

Ol 

5 

1 

'W 

i  i 

n  ,  exD 
ooeh  1 

L2 

+  ^ 

z 

,  ^  J 

j 

(4.16) 


The  last  factor  in  Eq.  (4.16)  constitutes  a  correction  factor  for  the  electron 
density  in  Eqs.  (4.1)  -  (4.14),  at  locations  close  enough  to  that  of  the  barrier 
(Fig.  4.1b).  Since  the  direction  of  the  y-coordinate  in  Eqs.  (4.1)  -  (4.13)  is 
more-or-less  aligned  with  that  of  the  barrier,  we  shall  use  this  correction  factor 
as  a  uniform  multiplicative  constant  to  be  applied  to  the  ambient  density  n^^^ 
wherever  appropriate  in  applications  (Sections  4.4  and  4.5). 


The  accuracy  of  this  procedure  is  not  immediately  clear.  As  one  moves 
outward  from  the  barrier,  $  becomes  less  negative  than  and  the  emitted- 
electron  contribution,  which  is  the  dominant  contribution  to  ng  near  the  barrier 
in  Eq.  (4. 15),  decreases  rapidly.  As  a  result,  it  becomes  less  clear  than  before 
whether  the  approximation  \d2<p/dx2 1  <<|32cf>/3y2|,  which  underlies  the  analysis 
presented  in  Eqs.  (4.1)  -  (4.13),  is  valid.  If  it  is  not,  then  |320/3y2|,  and 
therefore  |0|,  are  overestimated  by  this  procedure.  One  should  therefore  regard 
the  last  factor  in  Eq.  (4.16)  as  an  upper  bound  on  the  true  correction  to  <p 
resulting  from  emitt.ed-electron  space  charge.  In  representing  the  emitted 


electrons  as  Maxwellian,  we  have  lpnored  the  fact  that  the  velocity  distribution 

'  O  j 

of  backscattered  electrons  is  more  poorly  approximated  by  a  Maxwellian  and  also 
has  a  much  larger  average  energy  than  those  of  either  secondary  electrons  or 
photoelectrons  (Sternglass,  1954.  Fig.  8).  However,  at  the  incident-electron 
energies  of  Importance  In  our  applications  (Sections  4.4  and  4.5).  backscattered 
electrons  make  a  relatively  minor  contribution  to  the  total  emitted  flux. 

In  order  to  estimate  'I'g  using  the  results  presented  in  Eos.  (4.8)  -  (4.16)  and 
Figs.  4.2  and  4.3,  we  now  need  to  develop  a  way  of  estimating  the  ratio 
nwi/n  h  in  ^-9S-  (4*13)  and  (4.14).  Calculation  of  a  realistic  distribution  for 
ion  density  in  a  spacect  al  t  Wck.*,  in  the  presence  of  ion-drift  and  space-charge 
effects,  is  a  task  of  great  difficulty  (Sec.  4.4;  Sec  5).  A  very  simple 
approximation  is  the  "neutral  approximation"  in  which  electric-field  effects  on 
particle  motion  are  ignored.  For  a  plasma  having  ion  speed  ratio  S.  = 

U/(2kT./m.)  ,  this  approximation  readily  yields  the  result: 

-Sj  cos  0 

r^/n.^  =  7 r  ^  j  e  U  du  =  ^  1  -  erf  (S.  cos  0)  (4.17) 

-oo 

[see,  for  example,  Parks  and  Katx,  1983,  F.q.  (27)],  for  ion  density  on  a  convex 
surface,  where  0  is  the  angle  between  the  surface  normal  and  the  downstream 
direction.  Because  this  result  ignores  the  curvature  of  ion  orbits  caused  by 
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electric  fields,  it  grossly  underestimates  ion  density  on  (or  near)  downstream¬ 
facing  surfaces. 


To  obtain  a  better  estimate,  we  use  a  method  similar  but  not  identical  to  that 
of  Parks  and  Katz  (1983,  Sec.  4 ;  1985).  Here,  we  specifically  assume  that  the 
spacecraft  is  an  infinite  circular  cylinder  whose  axis  is  perpendicular  to  the  ion 
drift  direction.  We  assume  that  ions  arriving  at  the  downstream  axial  surface 
point  (r=r  ,  0  =  0  in  Fig.  4.1)  have  moved  through  a  potential  distribution  which 
has  the  form: 


0(r,0)  =  g(0)/r2 


(4.18) 


where: 


g(0)  = 


r 

s 


0o  (1  -  402/tt2),  for  1 01  <  rr/2 


^  0  otherwise 


(4.19) 


and  0O  is  the  potential  at  (r  =  r  ,  9  =  0).  As  in  Eq.  (4.3),  we  have  0O  <  0.  This 
choice  for  the  potential  is  different  than  that  of  Parks  and  Katz  (1983,  Eq.  40; 
1985,  Eq.  14).  It  is  a  special  case  of  the  form  0(r,0)  =  f(r)  +  g(0)/r2,  and  for 
potentials  of  this  form,  the  quantity: 

C  =  L3  +  mq  g(&)  (4.20) 
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is  a  constant  of  particle  motion  (Goldstein,  1980,  pp.  454-457).  Here,  m  and 

q  are  ion  mass  and  charge,  and  L  =  mrvg  is  ion  angular  momentum  in  the  (r,0) 

plane.  For  the  potential  given  by  Eqs.  (4.18)  and  (4.19),  this  feature  permits 

ion  orbits  arriving  at  (r  =  r  ,  0  =  0)  to  be  integrated  analytically.  The  density 

at  this  point  is  then  given  by  a  velocity-space  integral  whose  integrand  is  known 

analytically.  The  integral  itself  must  be  evaluated  numerically.  Details  are 

given  in  Appendix  C.  The  computer  program  used  for  doing  this  calculation  is 

listed  in  Appendix  D.  In  contrast  with  the  density  integration  of  Parks  and  Katz 

(1983,  1985),  ours  does  not  involve  an  inequality  (Parks  and  Katz,  1983,  Eqs. 

28-31,  cr  1985,  Eqs.  2-5).  The  potential  given  by  (4.18)  and  (4.19)  is  not  the 

same  as  the  one  given  by  Eqs.  (4.1  -4.14),  but  this  is  net  surprising  since  our 

treatment  is  not  self-consistent.  As  we  point  out  in  Sections  4.4  and  5,  a  self- 

consistent  treatment  of  our  problem  is  a  task  of  great  difficulty,  and  our 

approximate  treatment  serves  the  purpose  of  providing  a  simple  estimate  of 

wake  potentials.  In  the  potential  given  by  Eqs.  (4. 18-4. 19),  some  ion  orbits 

have  angular  turning-points  (Appendix  C  and  Fig.  4.4).  Figure  4.5  shows 

normalized  ion  densities  n.  =  n./n.  at  the  downs  -'cam  surface  point  (  r  =  r  , 

0  =  9),  in  the  presence  of  this  potential,  calculated  as  described  in  Appendices 

C  and  D.  these  values  can  be  used  to  estimate  the  ratio  n  ./n  ,  in  Eqs.  (4.131 

wi  ooen 

and  (4.14). 
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For  completeness,  we  also  develop  a  rough  correction  of  the  normalized  wake 
half-width  f  =  y  /A^  ,  for  effects  of  finite  speed  ratio.  We  do  this  for 
the  two  geometries  shown  in  Fig.  4.6.  In  terms  of  these  geometries,  we  do  this 
by  using  the  distance  BC,  rather  than  the  spacecraft  radius  or  half-width  OA  =  r  , 
for  Ymax-  We  assume  that  the  waKe  edge  AD  is  straight  l.see  Sec.  4.4)  and  has  a 
slope  of  -  1/S..  The  geometries  of  Figs.  4.6(a)  and  (b)  then  imply  the  results: 
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for  a  cylinder  or  sphere,  and: 
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(4.22) 


for  a  thin  disk  or  plate  normal  to  the  ion  drift  direction.  This  derivation  implies 
that  in  Eqs.  (4.1)  -  (4.13)  only,  the  y-coordinate  is  most  usefully  regarded  as 
distance  from  B  along  the  line  BC  in  Fig.  4.6(a)  or  (b).  As  long  as  at 
T.  ~  0.1  eV  is  the  dominant  ionic  constituent,  S.  is  large  enough  that  the 
corrections  given  by  (4.21)  or  (4.22)  are  unimportant  (Sec.  4.4). 


This  completes  the  development  of  our  method  for  estimating  the  potential  <p0 
at  the  downsteam  point  (r,6)  -  (r  ,0).  In  Sections  4.4  and  4.5,  we  apply  this 
method  to  wake-induced-barrier-effect  charging  of  the  Shuttle  Orbiter  and  of 
hypothetical  future  spacecraft  of  larger-  size. 
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4.4.  APPLICATIONS  TO  [  .C  W-PO:  .AR-ORRIT  SPACECRAFT  CHARGING; 
DISCUSSION 

We  examine  first  the  charging  of  the  Shuttle  Orbiter  in  the  auroral  plasma. 
Gur  primary  goal  is  to  find  out  if  there  exist  auroral-electron  energy 
distributions  which  are  not  energetic  enough  to  cause  high-voltage  charging  by 
overcoming  secondary-electron  and  backscattered-electron  emission  by 
spacecraft  surfaces,  but  which  can  cause  it  via  the  wake-induced-barrier 
mechanism.  The  range  of  auroral -electron  equivalent  temperatures  implied  by 
these  criteria  is  from  a  few  hundred  eV  to  a  few  keV,  with  the  upper  end  of  this 
range  being  strongly  dependent  on  the  secondary  and  backscattered  emission 
properties  of  the  surface  material  involved  (Sec.  4.1).  Within  this  range,  the 
most  important  remaining  criterion  is  a  small  enough  auroral-electron  Debye 
length,  which  implies  a  large  enough  auroral-electron  density. 

The  best  data  available  on  energy  distributions  of  auroral  electrons  are  those 
of  Gussenhoven  et  al  (1985)  and  Yeh  and  Gussenhoven  (1987),  obtained  using  the 
DMSP  F6  and  F7  satellites.  Gussenhoven  et  al  (1985.  Fig.  1 1)  present 
Maxwellian  fits  to  the  higher-energy  (above  a  few  keV)  portions  of  three  electron 
energy  distributions  observed  during  a  high-voltage  charging  event  on  DMSP  F7 . 
These  fits  yield  n_^  =  3.9,  3.2,  and  4.9  cm  ^  and  kTg^  =  10.1,  14.4,  and  4.2 
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keV,  resulting  in  Xq6^  values  of  378,  499,  and  218  m,  respectively.  For 
secondary  electrons  having  emission  temperature  =  3  eV,  applying  the 

density-correction  factor  A  +  A  -Ten/Ts^  from  EcF  (4-16)  gives  decreased  XQeh 
values,  respectively,  of  70,  84,  and  SO  m.  The  appropriate  value  of  r^  for  the 

Shuttle  Orbiter  depends  on  its  orientation.  Its  fuselage  radius  is  roughly  3m. 

Its  largest  wake  half-width  will  be  produced  when  its  wing  plane  is 
perpendicular  to  its  motion.  The  largest  circle  which  can  be  inscribed  inside 
its  cross-sect icr.  in  this  plane  has  a  radius  of  7.5m.  The  Orbiter’s  projected 
dielectric  area  in  the  same  plane  is  equivalent  to  that  of  a  circular  disk  of 
radius  11.5  m.  All  of  these  dimensions  are  much  smaller  than  the  Debye-length 
values  listed  above.  Equation  (4.9)  or  the  uppermost  curve  in  Fig.  4.3  now 
indicates  that  the  wake-induced -barrier  mechanism  will  induce  potentials  of  only 
a  few  volts  negative,  or,  with  the  density-corrected  Xq^  values  noted  above,  a 
few  tens  of  volts  negative.  The  largest  such  value  is  <p0  =  -  111  V, 
corresponding  to  kT0h  =  4.2  keV,  XDeh  =  50  m,  and  rs  =  1 1.5  m.  Evidently, 
smaller  values  of  Xpe^,  corresponding  to  greater  electron  densities  (Eq.  4.12), 
are  required  for  a  reliable  indication  of  high-voltage  charging  by  this 
mechanism. 

t eh  and  Gussenhover.  (1987,  Tables  2a-c  and  Figs.  1-2)  present  a  larger 
variety  of  auroral-zone  electron-energy  distributions.  Among  these,  two  result 
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ip.  Debye-iength  values  much  smaller  than  the  others.  These  are  the  "Type-i 

spectra"  given  in  their  Table  2a.  They  give  two  such  spectra.  These  are 

averages  of  spectra  which  produced  negative  DMSP  spacecraft  potentials  -<p 

greater  than  10CV  and  ir.  the  range  from  30V  to  70V,  respectively.  Each  of 

these  is  represented  by  two  different  Maxwellian  fits  below  and  above  a 

transition  energy  of  12.16  and  1  1.34  keV,  respectively.  For  the  first  spectrum, 

-  3 

the  parameters  characterizing  the  two  Maxwellian  fits  are  n  ,  =  58.30  cm 

o  el’ 

kT^,  =1.31  keV,  n^2  =  1-23  cm  ^  and  kT^  -  5.49  keV.  For  the  second, 
these  parameters  are  ngi  =  46.88  cm  ,  kTgi  =  1.14  keV,  neo  =  0.08  cm  , 

and  <Tg2  =  5.68  ke'v'.  The  most  important  difference  between  these  two  spectra 
:s  that  the  density  of  the  hotter  Maxwellian  is  much  larger  in  the  first  one.  As 

the  authors  point  out,  this  is  consistent  with  the  idea  that  the  higher-energy 
electrons  produce  combined  secondary  and  backscattered  yields  less  than  unity 
and  are  therefore  responsible  for  the  high-voltage  charging.  The  Debye  lengths 
^Deh  ^cr  “^ese  tw0  r>Fectra  are  35.2  m  and  36.7  m,  respectively,  based  on  the 
lower-energy  Maxwellian  fit;  the  presence  of  the  higher-energy  Maxwellian  has 
negligible  effect  on  these  values.  These  values  are  close  to  the  smallest  hct- 
electron  Debye  lengths  which  are  likely  to  be  encountered  in  the  auroral  zones 
(M.S.  Gussenhoven,  private  communication,  1987).  The  analysis  which  follows 
is  therefore  based  on  the  lower  energy  component  of  the  first  spectrum. 


F°r  XQeh  =  35.2  m  and  the  approximate  Shuttle  Orbiter  wake  half  widths  of 
7.5  m  and  11.5  m,  mentioned  above,  we  obtain,  using  either  F.q.  (4.10)  or  the 
uppermost  curve  in  Fig.  4.3,  the  surface  potentials  <p0  —  -29.7  V  and  -70.3V. 
Again  using  Eq.  (4.9)  with  the  density-correction  factor  from  Eq.  (4.16)  and  Tg 
=  3eV,  we  obtain  -  10.7m  and  the  surface  potentials  <p0  =  -273V  and 

-544V. 

The  variation  among  these  four  surface-potential  values  is  rather  large,  but 
their  peometric  mean  is  -133V,  so  if  we  take  this  as  the  most  realistic 
prediction  obtainable  using  our  analysis,  we  have  a  result  which  fulfils  the 
definition  of  "high-voltage"  charging  as  given  in  Sec.  1.  It  therefore  appears  that 
the  Shuttle  Orbiter  is  large  enough  to  undergo  high-voltage  charging  caused  by  the 
wa' ^-induced-barrier  mechanism,  at  least  for  some  combinations  of  spacecraft 
orientations  and  known  environmental  conditions.  For  the  same  conditions,  the 
same  methods  predict  that  in  order  to  achieve  (p0  =  -  1  00V  without  invoking  the 
density  correction  from  Eq.  (4.16),  a  spacecraft  half-width  of  14m  is  required, 
and  this  is  only  moderately  larger  than  the  Shuttle  Orbiter.  For  a  spacecraft 
half-width  of  60m,  the  corresponding  calculation  gives  <pc  =  -1000V;  however, 
this  latter  result  is  affected  strongly  by  ion-density  effects,  and  we  discuss  these 
later  in  this  Section. 
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One  consequence  of  this  charging  is  particularly  noteworthy.  Suppose  that 
the  downstream -facing  surfaces  of  a  spacecraft  are  mostly  dielectric,  but 
somewhere  near  the  center  of  the  downstream  surface  region  is  a  small  area  of 
exposed  conductor  which  is  grounded  to  the  spacecraft  frame,  and  this  in  turn  is 
grounded  to  exposed  conductive  surfaces  which  are  not  in  the  wake  and  therefore 
are  exposed  to  the  full  ram  flux  of  upstream  ions,  so  that  these  surfaces  will 
float  close  to  space  potential.  Evidently,  the  region  around  the  edges  of  the 
exposed  downstream  conductor  will  then  be  subjected  to  high  electrical  stresses, 
producing  a  prime  location  for  unwanted  electrical  discharges.  A  conductive 
surface  on  an  experiment  package  in  the  Orbiter’s  cargo  bay  is  an  example  of 
such  a  situation. 

In  the  above  derivation  of  spacecraft  potentials,  our  most  important 
approximation  was  the  neglect  of  the  ion  space-charge  density  in  the  wake.  The 
derivation  in  Sec.  4.3  and  Appendix  C  gives  us  the  means  to  estimate  this 
quantity.  In  order  to  do  this,  we  assume  that  the  ambient  plasma  has  a 
maximum  total  ion  or  electron  density  n  between  105/cm3  and  5  x  10s/cm3 
in  the  Shuttle  altitude  range  (Narcisi,  1973;  Samir  et  al,  1981,  Fig.  7; 

Swider,  1985,  rigs.  21-90  and  21-91).  The  data  presented  by  Rich  (1983, 
Fig.  3)  indicate  a  wider  range  of  possible  values.  We  assume  that  the  cold-ion 
const ituent  has  a  temperature  of  0. 1  eV.  We  further  assume  that  these  ions  are 
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either  O  or  H+  and  that  a  1  keV  hot-ion  constituent  may  or  may  not  be  present. 
For  a  spacecraft  orbital  speed  of  8  km/sec,  and  for  0+  at  0.1  eV,  H+  at  0. 1 
eV,  O  at  1  keV,  and  H  at  1  keV,  the  ion  speed  ratio  S.  has  values  of  7.31, 
1.83,  0.0731,  and  0.0183,  respectively;  the  last  two  of  these  are  effectively 
zero.  In  calculating  the  above-mentioned  spacecraft  potentials,  we  used  a  hot- 

electron  amoient  density  value  n  ^  of  58.3  cm  .  This  is  smaller  than  the 

above-mentioned  total  n  values  by  factor's  of  1.7  x  103  and  8.G  x  103, 

ooe  ' 

respectively.  In  order  that  our  predicted  spacecraft  potentials  not  be  greatlv 
decreased  by  effects  of  ion  space-charge,  the  ion  density  n^.  in  the  wake  must 
be  reduced  below  ambient  by  an  even  larger  factor,  say  10s.  We  now  see  from 
Fig.  4.5  that  for  S,  =  7.31,  such  a  reduction  requires  -y0  =  -e<p0/WT^  <  about 
1000,  i.e.,  -  ct>o  4  about  100  V.  For  larger  values  of  -(£0,  ion-density  effects 
will  begin  to  alter  the  value  of  cp0  significantly.  Using  the  same  n^  values 
together  with  Eq.  (4.14]  and  Fig.  4.5,  we  obtain  limiting  <fi0  values  of  -311  V 
and  -198  V,  respectively,  in  the  limit  of  large  spacecraft  size.  If  we  again  use 
the  density-correction  factor  from  Eq.  (4.16)  with  Ts  =  3  eV,  these  4>0  values 
increase  to  -727  V  and  -395  V,  respectively. 

If  any  substantial  fraction  of  H+  is  present,  or  any  substantial  ion  fraction 
having  T;  ~  1  keV,  then  Figs.  4. 2-4. 5  indicate  that  these  ions  will  suppress 
wake- induced-barrier-effect  charging. 
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At  300  km  altitude,  measured  midlatitude  H+  concentration  values  of  64/cm3 
and  2.0  x  103/cm3  are  given  by  Narcisi  (1973,  Figs.  5  and  6);  see  also  Swider 

(1985,  Figs.  21-90  and  21-91).  Two  measurements  of  2.8  x  l02/em3  at  300 
km  have  been  given  by  Samir  et  al  (1981,  Figs.  7a  and  7b).  A  measurement 
close  to  this  value  has  also  been  given  by  Rich  (1983,  Fig.  17).  For  the  above- 
mentioned  T.  and  S.  values  of  0.1  eV  and  1.83,  Fig.  4.5  indicates  that  the  wake- 

point  ion  density  at  <po  =  ~  100V  is  reduced  below  the  ambient  value  by  a  factor 

-2  +  2  2 
of  4.3  x  10  .  This  gives  wake-point  H  concentrations  of  2.8/cm'3,  86/cm  , 

3 

and  12/cm  ,  respectively.  The  first  and  third  of  these  are  substantially  below 

-3 

the  n  ^  value  of  58.3  cm  mentioned  above,  indicating  that  in  these  cases,  not 

enough  H+  is  present  to  suppress  wake-induced  charging.  Furthermore,  the  same 
measurements  show  smaller  concentrations  at  lower  altitudes.  Narcisi 

(1973)  also  presents  data  showing  decreases  in  both  0+  and  H+  concentrations  at 
auroral  latitudes  compared  to  midlatitudes.  In  addition,  one  does  not  necessarily 
find  kilovolt  ions  when  kilovolt  electrons  are  present  (M.S.  Gussenhoven,  private 
communication,  1987).  Furthermore,  the  wake-point  ion-density  values 
displayed  in  Fig.  4.5  may  be  too  large,  for  reasons  discussed  following  Eq. 

(Cl  7).  For  all  these  reasons,  it  therefore  appears  that  the  wake-induced-barrier 
mechanism  can  produce  high-voltage  charging  on  spacecraft  at  least  as  large  as 
the  Shuttle  Orbiter. 
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The  foregoing  predictions  involve  considerable  uncertainty.  The  most 

important  causes  of  this  are:  (1)  the  lack  of  simultaneity  of  the  above-mentioned 

measurements  of  n  ,  ,  n  ,  and  n  ,  (2)  the  crudity  of  our  correction  for 

coen  H+  Q+ 

emitted-electron  density,  which  involved  the  use  of  the  last  factor  in  Eq.  (4.16) 
as  a  spatial ly-uniform  factor,  and  (3)  the  crudity  of  our  method  for  calculating 
the  nwi  values  presented  in  Fig.  4.5.  The  major  source  of  crudity  in  item  (3) 
is  our  assumption  of  a  model  potential  [Eqs.  (4.18)  and  (4.19)]  which  probably 
approximates  the  actual  self-consistent  potential  distribution  badly  enough  to 
produce  important  errors.  To  do  better  is  a  task  of  great  computational 
difficulty.  To  calculate  downstream -point  surface  potentials  exactly  requires 
knowing  not  just  one  value  of  ion  density  but  rather  the  entire  ion-density 
distribution  around  the  spacecraft,  and  this  is  coupled  with  the  entire  potential 
distribution  in  a  strongly  nonlocal  and  nonlinear  way  through  Poisson’s  equation 
and  the  particle  dynamics.  If  most  of  the  ions  are  0+  at  ~  0.1  eV,  the  ambient 
ion  velocity  distribution  is  very  anisotropic.  The  ram  energy  of  an  0+  ion  at  8 


km /sec  is  about  4  eV,  and  ions  entering  a  high-voltage  wake  are  accelerated  to 
kinetic  energies  much  larger  than  this.  The  ion  velocity  distribution  in  such  a 
wake  will  therefore  be  very  "beam-like"  and  the  resulting  ion-density  distribution 
is  likely  to  be  very  sensitive  to  small  changes  in  spacecraft  geometry  and 
external  conditions,  and  also  to  small  changes  in  the  potential  distribution  of  the 
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kind  encountered  in  an  iterative  numerical  calculation.  The  latter  type  of 
sensitivity  may  destabilize  such  a  calculation.  A  numerical  calculation  by 
Laframbcise  and  Parker  (1987),  of  the  surface-potential  distribution  on  a 
nonemitting  dielectric  cylinder  in  a  collisionless  plasma  crossflow  containing  a 
hot-electron  constituent,  illustrates  many  of  these  features.  A  description  of 
their  work  also  appears  in  Sec.  5.  A  more  thorough  investigation  of  the  wake- 
induced-barrier  effect  will  require  numerical  calculations  in  which  all  of  these 
difficulties  will  be  encountered. 
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4.5.  APPLICATIONS  TO  BEAM-INDUCED  SHU  1  TLE  ORBITER  CHARGING; 
DISCUSSION 

Our  second  investigation  concerns  charging  of  the  Shuttle  Orbiter  during 
electron-beam  emission  by  the  SEPAC  Spacelab-1  experiment.  Sasaki  et  al 
(1986)  report  that  charging  of  the  Orbiter  due  to  beam  emission  was  strongly 
dependent  on  the  orientation,  of  the  Orbiter  with  respect  to  its  velocity  vector. 
For  beam  emission  currents  ^  100  mA,  the  Orbiter’s  frame  potential  was 
observed  (see  their  Fig.  10)  to  reach  the  beam  acceleration  voltage  when  its 
cargo  bay  faced  downstream,  but  not  otherwise.  Herein,  we  propose  an 
explanation  for  this  observation,  based  on  the  wake-induced-barrier  mechanism. 
The  basic  premise  is  that  if  the  "upper"  side  of  the  Orbiter  faces  downstream, 
then  a  iarge  region  of  space  surrounding  the  cargo  bay  will  be  depleted  of  ions. 

If  the  wake- induced-barrier  mechanism  then  produces  negative  potentials  on  the 
Orbiter’s  upper  wing  surfaces  and  other  wake-region  dielectric  surfaces,  the 
region  of  positive  potentials  (the  "potential  well  for  electrons"),  which  must 
surround  the  exposed  conductive  surfaces  associated  with  the  SEPAC  experiment 
in  the  cargo  bay,  will  then  be  very  constricted  in  spatial  extent,  and  the  process 
of  electron  return-current  collection  will  be  partly  inhibited,  requiring  larger 
positive  potentials  on  these  surfaces  in  order  to  collect  the  electron  current 
required  for  current  balance. 
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In  order  to  investigate  this  situation  in  detail,  we  use  data  presented  by 
Taylor  et  al  (1985).  Plate  3  of  their  paper  shove  the  response  of  the  SEPAC 
energetic-electron  analyzer  during  a  time  interval  which  included  six  pulses  from 
the  SEPAC  electron  gun.  The  duration  of  each  pulse  was  4  seconds.  This  time 
interval  formed  part  of  the  SEPAC  FO-2  experiment  sequence,  during  which  the 
Orbiter’s  cargo  bay  faced  downstream.  We  analyze  the  electron  spectrum 
measured  during  the  second  pulse,  during  which  the  electron  beam  current'  and 
emission  energy  were  70  mA  and  5  keV.  During  this  pulse,  their  Plate  3  shows 
a  sharp  peak  near  1  kV  in  the  spectrum  measured  by  the  energetic-electron 
analyzer,  indicating  Orbiter  frame  charging  to  that  voltage.  The  same  spectrum 
also  shows,  during  all  six  pulses,  an  enhanced  population  of  electrons  up  to  the 
upper  end  of  the  detector  range  at  15  kV,  i.e.  well  above  beam  energy.  These 
electrons  presumably  were  energized  by  a  wave-particle  interaction  caused  by  the 
presence  of  the  beam  (Taylor  et  al,  1985).  The  population  above  1  kV  decreases 
monotonically  with  increasing  energy. 

We  need  to  find  out  whether  these  energized  electrons  can  produce  a 
significant  voltage  on  downstream -facing  dielectric  surfaces  via  the  wake- 
induced-barrier  mechanism.  To  do  this,  we  must  estimate  their  "ambient" 
density,  i.e.  that  at  space  potential.  We  do  the  latter  by  fitting  a  Maxwellian 
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distribution  to  the  above-mentioned  spectrum.  This  spectrum  shows  energy- 
differential  fluxes  $  =  6.31  x  10s  and  1.585  x  10s  electrons/cm2  sec  steradian 
eV  at  energies  E  of  about  3.59  and  12.2  keV  with  respect  to  spacecraft  frame 
potential,  respectively.  Assuming  that  these  electrons  have  an  isotropic  velocity 
distribution,  we  have: 
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(Prokopenko  and  Laframboise,  1977,  1980),  where  dJeo/dE  is  the  energy- 
differential  flux  onto  one  side  of  an  arbitarily-oriented  surface  element,  and  f  is 
the  electron  velocity  distribution  function.  By  fitting  a  Maxwellian  distribution 
of  the  form: 
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to  the  resulting  values  of  fg,  we  obtain  ng  =  18.6  electrons/cm3  and  kTg  = 

3.32  keV.  This  is  the  value  of  ng  at  spacecraft  frame  potential,  so  it  must  be 
multiplied  by  exp  (-1.0/3.32)  to  obtain  ng  at  space  potential.  Doing  this  yields 
the  "ambient"  value  ng  =  13.8  electrons/cm3.  Again  using  F.q.  (4.10)  and  our 
estimates  of  7.5  and  1 1.5  m  for  Orbiter  wake  half-width,  together  with  these 
values  of  ng  and  kTg,  we  obtain  <p0  =  -7.0  or  -16.6V,  respectively,  as  estimates 
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for  the  wake-induced-barrier  potential  for  this  situation.  Again  using  the 
electron-density  correction  factor  from  Eq.  (4.16)  with  kTs  =  3  eV  as  in  Sec. 
4.4,  we  obtain  <p0  =  -117  or  -266  V.  As  in  Sec.  4.4,  we  propose  the  geometric 
mean  of  these  four  values,  -43.6  V,  as  a  rough  "best  guess"  for  the  actual 
potential  achieved.  As  in  Sec.  4.4,  these  four  estimates  cover  a  large  range  of 
potentials.  The  smallest  among  them  is  not  much  larger  than  the  potential 
difference  of  a  few  volts  expected  between  the  barrier  and  the  surface  potentials 
(Sec.  4.2),  but  they  still  serve  to  indicate  that  when  the  Orbiter’s  cargo  bay 
faced  downstream  during  the  SEPAC  experiment,  the  wake-induced-barrier 
mechanism  probably  caused  at  least  marginally  negative  potentials  to  occur  over 
most  of  the  Orbiter’s  downsteam-facing  dielectric  surfaces.  The  resulting 
steepening  of  the  "potential  well  for  electrons"  around  the  SEPAC  experiment 
would  then  restrict  electron  collection,  by  causing  a  breakdown  of  "orbit- 
limitation"  of  this  collection  (Mott-Smith  and  Langmuir,  1926;  Laframboise  and 
Parker,  1973)  and  this  would  force  the  spacecraft  frame  potential  to  become 
more  positive  in  order  to  collect  the  r>  quired  return  current. 

In  contrast  with  this,  one  can  readily  infer  that  when  the  Orbiter’s  cargo  bay 
and  upper  surfaces  faced  upstream  rather  than  downstream,  potentials  on  a 
larger  portion  of  these  surfaces  would  have  been  positive  with  respect  to  space, 
by  amounts  which  increased  with  decreasing  distance  from  the  SEPAC  electron 
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gun.  To  show  this,  we  proceed  as  foliows.  Sasaki  et  al  (1986)  assume  that  the 
ambient  plasma  surrounding  the  Orbiter  had  a  temperature  of  1000  K  and  a 
density  of  iOs/cm3.  The  corresponding  random  electron  flux  is  0.79  mA/m2. 

In  comparison,  the  random  flux  for  the  Maxwellian  parameters  calculated  just 
above  is  0.02  mA/m2,  and  this  is  very  much  smaller.  In  cases  where  the  cargo 
bay  faced  upstream,  densities  of  ambient  ions  near  it  would  be  much  larger,  and 
the  wake-induced  potential  barrier  would  be  absent.  Using  the  backscattered- 
electron  and  secondary-electron  yield  data  for  fused  silica,  the  Shuttle  dielectric 
surface  material,  given  by  Prokopenko  and  Laframboise  [1980,  Eq.  (9)], 
Laframboise  et  al  (1982,  Table  1),  and  Yeh  and  Gussenhoven  (1987,  Fig.  3, 
rescaled  for  fused  silica),  we  obtain  that  the  total  emitted-eicctron  flux  exceeds 
the  incident-electron  flux  for  incident  -electron  energies  from  40  eV  to  23.3  keV. 
Therefore,  the  impingement  of  monoenergetic  electrons  anywhere  in  this  energy 
range  constitutes  a  net  positive  rather  than  negative  current,  and  if  no  other 
current  sources  were  present,  this  current  would  drive  a  surface  exposed  to  it  to 
a  positive  potential  large  enough  that  the  beam  would  impact  at  a  kinetic  energy 
of  23.3  keV;  this  energy  is  called  the  "second  crossing"  of  the  emission  yield 
curve  for  this  material  (Leung  et  al,  198  1L 

For  a  beam  emitted  at  5kV  with  respect  to  spacecraft  ground,  this  implies 
that  dielectric  surfaces  exposed  to  it  would  float  at  18.3  kV  with  respect  to 
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spacecraft  ground,  assuming  that:  (a)  it  remained  nearly-enough  monoenergetic, 

[b)  its  impacting  current  density  or,  surfaces  close  to  space  potential  were  enough 
to  overcome  the  negative  current  density  from  the  low-energy  ambient  electrons, 
and  (c)  the  emitted  electrons  all  or  almost  all  escaped.  If  these  dielectric 
surfaces  thereby  became  the  most  positively-charged  locations  on  the  spacecraft, 
assumption  (c)  would  be  violated,  since  err- *t  fed  cle-d  .vr.©  generally  tend  to 
migrate  to  the  most  positive  electrode  available,  and  this  would  prevent 
dielectric  surfaces  from  acquiring  voltages  much  more  positive  than  those  of  the 
exposed  conductors.  These  would  be  at  spacecraft  frame  potential,  which  would 
be  between  0  and  5  kV  positive  with  respect  to  space. 

To  find  out  if  returning  beam  fluxes  can  ever  overcome  the  cold-electron 
ambient  flux,  we  note  that  the  largest  electron  current,  emitted  by  the  SEP  AC 
gun  was  300  mA.  If  all  of  this  current  were  to  return  and  impact  uniformly 
over  an  area  equal  to  the  Orbiter’s  projected  wing-plane  cross-section  of  418 
m2,  the  resulting  electron  flux  would  then  be  0.72  rr,A/m2,  and  this  is  close  to 
the  cold-electron  random  flux  value  of  0.79  mA/m2  noted  above.  Because 
ambient-electron  current  is  also  present,  not  all  of  the  beam  current  will 
return,  even  when  the  spacecraft  frame  is  charged  to  beam  voltage,  but  the 
returning  beam  current  is  likely  to  be  more  concentrated  at  smaller  distances 
from  the  SEPAC  experiment,  so  this  comparison  suggests  that  a  relatively 
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extended  neighborhood  of  the  SEPAC  experiment  existed  within  which  dielectric 
surfaces  charged  to  positive  voltages  close  to  the  beam  emission  voltage,  at 
times  when  the  cargo  bay  was  exposed  to  ambient  upstream  ions.  The  resulting 
enlargement  of  the  positive-potential  region  surrounding  the  SEPAC  experiment 
should  then  have  permitted  easier  beam-current  return  to  the  spacecraft  frame, 
thereby  tending  to  prevent  the  spacecraft  frame  from  charging  to  the  beam 
emission  voltage,  and  this  appears  to  support  our  explanation,  proposed  herein, 
of  the  SEPAC  observations  of  spacecraft -frame  voltages.  The  edges  of  this 
region  might  be  expected  to  fluctuate  considerably,  and  this  may  explain  some  of 
the  time-dependence  reported  in  the  observations  of  Taylor  et  al  (1985)  and 
Sasaki  et  al  (1 986) . 

It  might  be  expected  also  that  when  the  cargo  bay  faced  upstream,  the 
resulting  increase  in  the  density  of  both  ambient  ions  and  electrons  near  it  would 
increase  the  tendency  of  the  beam  to  be  scattered  by  beam-plasma  instabilities 
near  the  Orbiter,  and  this  would  further  facilitate  beam-current  return. 

However,  no  such  tendency  is  apparent  when  one  compares  the  electron  spectra 
shown  in  Plate  3  (cargo  bay  facing  downstream)  and  Plate  T  (cargo  bay  facing 
upstream)  of  Taylor  et  al  (1985). 
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5.  PRELIMINARY  RESULTS  FROM  \  NUMERICAL  LOW-POLAR -ORBIT 
CHARGING  SIMULATION 
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Section  we  present  preliminary  results  from  numerical  simulation  work  whose 
initial  goal  wan  to  verify  this  suggest ior:.  This  work  uses  an  adaptation  of  a 
rrevmus! v-exislirg  Gmuiaticr.  program  written  by  one  of  us  -Parker,  19831  for  ar. 
wfirCe  cylindrical  soareoraft  geometry  in  a  collisionless  plasma  crossflow.  If  the 
spacecraft  surface  is  entirely  dielectric,  then  at.  progressively  increasing  values  cf 
R.  ore  might  expect  surfaces  in  the  spacecraft  wake  to  be  the  first  to  undergo  high  - 
vo ;  tag'-1  n’wrgirg,  fell  :wed  by  sideways-facing  and  finally  frontward  -facing  surfaces, 
a  :  shown.  srhemat.maLy  ir  Fig.  5.1.  The  ruin  erica!  results  show  that  at  least,  cne 
implied  feature  vf  Fig.  5.1,  namely  the  apparent  monotonic  progression  to  larger 
■~.egr : ve  voltage?  an  ire  moves  arouse  the  surface  ci  the  cylinder  toward  the  wake 
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study  included:  Debye  length  Ap/spacecraft  radius  r  -0.001,  ion  speed  ratio 

(  =  drift  speed  U/ion  most-probable  thermal  speed  v'2kT./mp  -  8,  ions  0+  at  a 
temperature  of  0.2  eV,  cold  electrons  at  0.2  eV,  and  hot  (auroral)  electrons  at  5 

keV.  We  have  again  made  the  assumptions  (1),  (2),  and  (3)  listed  in  Section  2. 

The  finite-element  discretization  used  for  position  space  is  indicated  in  Figs. 

5.2.  and  S.3.  "Inside-out"  ion  orbit-following  was  used  for  calculating  ion  surface- 
current  densities.  The  ion  velocity-space  discretization  used  was  similar  to  that 
described  by  Parker  (1977).  To  achieve  sufficient  accuracy  in  ion  flux  values  at 
wake-side  surface  points  required  use  of  1024  incident  ion  directions  at  each  of  32 
ion  energy  levels  at  each  such  point.  The  iteration  was  started  using  an  ion  density 
distribution  based  on  assuming  that  ions  behaved  as  neutrals,  and  iterating  until  the 
surface  potential,  surface  ion  current  density,  sheath  potential  and  sheath  electron 
density  distributions  all  converged.  Electron  densities  and  currents  were  described 
by  superpositions  of  Boltzmann  factors.  Each  such  calculation  took  about  6  hours  on 
the  AFGL  Cyber  850  computer.  The  next  intended  step  was  to  calculate  an  ion 
density  distribution  corresponding  to  the  resulting  sheath  potential  distribution,  then 
"freeze"  th.s  i on  density  distribution  and  iterate  the  other  quantities  involved  as 
before,  and  so  on.  To  carry  this  procedure  to  convergence  would  probably  take  20  to 
30  hours  on  the  same  computer,  so  we  have  not  so  far  made  such  calculations,  and 
the  results  presented  here  are  therefore  based  on  the  "neutral  approximation"  for  ion 
densities,  but  on  ion  orbit -following  for  ion  current  collection  on  surfaces. 

Preliminary  results  from  these  calculations  are  shown  in  Fig.  5.4,  for  hot- 
electron  to  total  ion  flux  ratios  R  [See  Lq.  (2.8)]  equal  to  0.1,  0.25,  0.5,  and 
5.0.  On  the  front  end  sides  of  the  cylinder  (0  4  80°),  surface  potentials  seem  to 
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increase  more-or  less  u ionotonicaiiy  as  one  moves  away  from  the  front  (0  =  180°), 
except  for  some  relatively  small  oscillations  which  are  probably  spurious.  One 
possible  cause  of  these  is  loo  coarse  a  position  space  discretization.  Successive 
modifications  of  the  computational  grid  to  reduce  such  errors  resulted  in  the  final 
grid  shown  in  Fig.  5.3,  in  which  nodes  are  placed  very  densely  just  outside  its 
surface,  and  the  density  of  nodes  is  also  increased  in  the  region  80°  <  G  <  100°. 
Other'  possible  causes  include  incomplete  convergence  or  insufficiently  fine  sampling 
of  the  ion  velocity  space.  Even  with  32,768  orbits  followed  per  surface  point,  very 
few  of  these  will  connect  back  to  the  ambient  plasma  with  an  ambient  velocity  close 
to  that  of  the  heavily -populated  part  of  the  ion  distribution.  This  strongly  suggests 
that  in  future  calculations,  outside-in  orbit-following  needs  to  be  used  for  defining  the 
most  important  region(s)  of  velocity  space  for  incoming  ion  orbits,  perhaps  combined 
with  inside-out  orbit  following  for  doing  detailed  sampling  of  these  regions. 

However,  on  rearward  surfaces  (0  s  80°),  there  is  some  non-monotonicity  which 
appears  to  be  real.  For  F<  =0.1,  0.25,  and  0.5,  there  are  three  surface-potential 
maxima,  one  at  the  rearmost  point  (0  =  0°)  and  one  on  each  side  of  it.  For  R  = 

0.  1  and  0.25,  these  features  were  almost  unchanged  (the  potentials  of  the  maxima 
changed  by  less  than  2%)  when  as  few  as  4,096  ion  orbits  (4  ion  energy  levels)  were 
used,  and  this  attests  to  their  reality,  f  or  R  =  5,  we  obtain  four  maxima,  two  on 
either  side  of  the  rearmost  point.  The  minima  between  these  oeaks  corresDond  to 

i  i  i 

the  ion  "hot  spots'",  or  deposition  poim  ;  of  highly  beam -like  ion  populations,  men¬ 
tioned  above.  Most  real  situations  would  not  possess  the  symmetry  about  the  fore-aft 
line  ww  d,  our'  problem  does,  but  our  results  nonetheless  indicate  that  high  voltage 
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This  is  related  also  to  the  fact  that  the  ambient  ion  distribution  is  a  highly-direc¬ 
ted  (S.  =  8)  one.  A  situation  with  ion  drift  is  very  different  than  one  without.  A 
drifting  distribution  is  not  an  equilibrium  one.  The  drift  provides  the  ions  with 
thermodynamic  free  energy  which  can  support  a  much  greater  variety  of  phenomena 
than  in  the  nondrifting  case,  including  self-excited  oscillations  (Krall  and 
Trivelpiece,  1973,  Chs.  2  and  9).  For  related  reasons,  the  drift  also  makes  possible 
ion  focusing  effects  including  those  discussed  here.  These  effects  depend  on  the  fact 
that  when  the  ambient  ion  velocity  distribution  is  anisotropic,  as  it  is  very  greatly 
when  S.  =  8,  a  knowledge  of  the  phase-space  boundary  separating  orbits  which  connect 
back  to  infinity  from  those  which  do  not,  no  longer  gives  complete  information  about 
the  ion  velocity  distribution  at  any  finite  location.  In  the  steady-state  treatment 
described  here,  self-excited  oscillations  cannot  be  treated.  However,  oscillations 
observed  in  the  outer  portions  of  the  Orbiter’s  wake  involve  density  fluctuations  of 
only  a  few  per  cent  (Murphy,  1985;  Murphy  et  al,  1986),  so  they  probably  do  not 
have  much  effect  on  the  steady-state  wake  properties  studied  here. 

Beam-like  ion  constituents  have  important  implications  for  more  complicated 
situations,  such  as  the  interior  of  the  Orbiter’s  car01  , .  As  our  work  has 
indicated,  beam  impact  points  will  be  very  sensitive  to  details  of  the  potential  distri¬ 
bution  on  or  near  a  spacecraft.  Ions  may  come  around  a  corner  of  the  spacecraft, 
just  miss  one  equipment  box,  and  impact  on  one  small  corner  of  the  next  box,  which 
may  then  come  to  a  potential  very  different  than  those  of  its  surroundings. 

In  addition,  more  than  one  ion  beam  may  be  present  at  some  points  in  the  sheath. 
Beamsmay  come  around  corners  from  opposite  sides  of  a  spacecraft,  and  beams  may 
also  be  present  which  have  circled  the  spacecraft  one  or  more  times.  Ion  velocity 
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distributions  as  a  functio: 


i  of  direction  may  therefore  contain  many  narrow  "peaks" 
and  "valleys'-.  Simulation  of  a  very  derailed  kind,  which  will  tax  the  capabilities  of 
the  fastest  computers  presently  available,  will  be  necessary  to  resolve  such 
situai ions. 

In  Fig.  2.2  we  have  also  plotted  the  largest  value  of  negative  surface  potential  for 
each  of  the  four  values  of  R  shown  in  Fig.  5.4,  for  comparison  with  our  simple 
prediction  from  Sec.  2.  These  values  are  shown  as  circled  points  in  this  Figure. 
These  results  qualitatively  resemble  the  * heoretical  curves  in  the  same  Figure,  and 
also  our  expected  behaviour  shown  in  Fig.  5.1,  in  that  they  appear  to  show  almost  no 
dependence  on  R  until  a  "threshold"  value  is  reached,  and  then  they  show  a  rapid 

increase.  Furthermore,  this  increase  appears  to  begin  at  a  substantially  lower  R 
value  than  for  these  curves,  as  one  expects  since  the  curves  are  for  whole-body 
charging  and  the  plotted  points  are  for  charging  of  electrically-isolated  surfaces  in 
the  wake. 


One  feature  of  these  results  is  unexpected:  the  "pre-threshold"  wake  charging, 
instead  of  being  close  to  zero,  is  already  several  hundred  volts.  However,  this 
apparent  threshold  is  not  the  true  one.  High-voltage"  charging  conventionally  means 
that  surface  potential  exceeds  100V  in  magnitude.  For  our  conditions,  and  for  a 
surface  potential  of  -100V,  the  hot-eieetron  flux  exceeds  the  cold-electron  flux  when 
R  >  1.2  x  !  O'21'’!  Also,  for  R  -  0.1  (the  left  margin  of  Fig.  2.2),  if  we  calculate 
down, stream -point  potential  using  the  neutral  approximation  for  ion  flux  [Tsien,  1946 
parks  arid  Kat7,  1983,  F.q.  4j,  we  uc'.ain  'pr,  -  -3.3  x  !0'V.  For  R  -  0  (cord 

3  orm!  a  su:  Far  calculation  yields  <p  -  -14.8V.  It  is  evident  that  all  the 
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shown  in  Fig.  2.2  are  ’post-threshold"  in  the  sense  that 


they  involve  a  current  balance  primarily  between  the  hot  electron  flux  and  wake- 
regior.  ion  fluxes  that  are  already  heavily  modified  by  mbit  curvature  in  strong  wake- 
region  electric  fields.  This  is  evidently  a  situation  in  which  even  a  small  amount  0f 
ambient  HT  can  be  expected  to  produce  a  large  decrease  in  the  magnitudes  of  wake- 
region  potentials.  In  Section  *4.4,  we  examined  H  effects  on  wake-region  potentials 
using  the  model -potential  calculation  of  ion  densities  developed  in  Section  4.3  and 
Appendices  C  and  D,  and  we  made  a  preliminary  determination  that  H+  densities  were 
too  small  in  most  low-polar-orbit  conditions  to  cause  a  significant  reduction  in  wake- 
mduced-barrier-effect  potentials.  In  the  situation  studied  in  this  Section, 
downstream  H+  fluxes  rather  than  densities  are  the  most  important  issue.  We  intend 
to  include  H+  effects  in  future  calculations. 
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6. 


CONCLUSIONS 


In  Section  1  of  this  Report,  we  reviewed  the  main  differences  between  the 
plasma  environments  in  geosynchronous  orbit  and  low  polar  orbit  with  respect  to 
nigh- voltage  charging  situations,  and  we  gave  reasons  why  the  task  of  predicting 
high-voltage  charging  or  related  environmental  effects  on  a  spacecraft  is  more 
complicated  in  the  low-polar-orbit  situation. 

In  Section  2,  we  developed  a  simple  rough  estimate  of  the  environmental 
conditions  required  for  overall  charging  of  a  large  spacecraft  in  low-orbit 
auroral-zone  conditions.  The  results  indicate  that  for  any  given  spacecraft, 
surface  potentials  are  likely  to  depend  more  strongly  on  the  ratio  of  ambient  flux 
of  high-energy  electrons  to  that  of  all  ions  than  on  any  other  environmental 
parameter.  This  prediction  has  been  corroborated  experimentally  by  results  of 
Gussenhoven  et  al  (1985)  obtained  using  the  DMSP  F6  and  F7  satellites;  see 
their  Fig.  7  and  associated  discussion,  and  also  Laframboise  and  Parker  (1986, 
1987). 

In  Section  3,  we  presented  an  approximate  theoretical  treatment  of 
suppression  by  magnetic  fields  of  the  escape  of  electrons  emitted  from 
negative! y-charged  spacecraft  surfaces.  Our  treatment  of  this  topic  is  intended 
to  serve  a  variety  of  purposes.  Its  results  provide  a  "parameterization11  of  the 
electron -escape  problem  which  oar,  oo  used  as  part  o:  large  spacecraft-charging 
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simulation  programs,  and  this  may  permit  important  reductions  in  the  cost  of 
using  such  programs.  Our  results  also  permit  simple  estimates  of  the 
importance  of  magnetic  electron-escape  suppression  in  particular  situations. 
However,  their  most  important  use  may  be  to  permit  gaining  an  understanding  of 
an  effect  which  appears  likely  to  dominate  the  high-voltage  charging  of  mostly- 
dielectric  large  spacecraft  such  as  the  Shuttle  Orbiter  in  low  polar  orbit.  The 
nature  of  our  results  strongly  implies  that  on  such  spacecraft,  local  charging, 
especially  on  downstream -facing  surfaces  nearly  parallel  to  B,  may  occur  in 
ionospheric  conditions  which  do  not  produce  overall  charging.  These  conditions 
may  include  lower  equivalent  auroral-electron  temperatures  and  also  daylit  parts 
of  the  orbit.  This  question  demands  experimental  investigations  and  detailed 
numerical  simulations. 

Our  results  in  Sections  3. 3-3. 5  also  show  "sensitivity"  effects  which  have 
disturbing  implications.  These  effects  imply  that  obtaining  realistic  predictions 
of  high-voltage  polar-orbit  chaining  in  geometrically  complicated  situations  is 
likely  to  prove  a  computational  task  of  great  difficulty.  An  experimental 
investigation  of  magnetic-field  effects  on  secondary-electron  and  photoelectron 
escape  from  negatively-charged  surfaces  would  be  of  great  value,  and  we  hope 
that  the  work  presented  in  Section  3  will  stimulate  such  efforts. 


In  Section  4,  we  developed  theory  which  permits  an  approximate  calculation 
of  potentials  produced  on  wake-region  spacecraft  surfaces  by  the  wake-induced- 
barrier  effect.  In  Section  4.4,  we  have  used  this  theory  to  estimate  surface 
potentials  on  large  spacecraft  covered  mostly  with  dielectric,  including  the 
Shuttle  Orbiter,  in  auroral-plasma  conditions.  In  Section  4.5,  we  have  proposed 
an  explanation,  involving  the  wake-induced-barrier  effect,  for  observations,  made 
during  electron-bearn  emissions  from  the  SEP  AC  Spacelab-1  experiment,  that  the 
Orbiter’s  frame  charged  to  beam  emission  voltage  when  the  emitter  faced 
downstream  but  not  upstream. 

In  auroral-plasma  conditions,  we  predict  that  the  largest  wake-induced- 
barrier-effect  potentials  will  occur  when  the  Orbiter’s  wing  plane  is 
perpendicular  to  its  orbital  motion  and  the  hot  (£  1  keV)  electron  constituent 
has  the  largest  possible  density  (Eq.  4.12).  For  the  largest  such  density  so  far 
observed,  we  predict  potentials  more  than  100V  negative  on  surfaces  located 
near  the  wake  center.  For  larger  spacecraft  in  the  same  plasma  conditions,  we 
predict  potentials  several  hundred  volts  negative. 

In  contrast  with  the  situation  for  low-polar-orbit  charging  already  observed  on 
the  DMSP  F6  and  F7  satellites  (Gussenhoven  et  al,  1985;  Yeh  and  Gussenhoven, 
198  7),  wake-induced-barrier-effect  charging  does  not  require  the  ambient  hot 
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electrons  to  be  energetic  enough  that  their  incident  flux  overcomes  secondary 
and  backscattered  emission  fluxes,  because  these  are  suppressed  by  the  potential 
barrier  involved  in  this  mechanism.  As  in  the  case  of  suppression  by  B  effects, 
this  extends,  to  lower  equivalent  temperatures,  the  range  of  auroral -electron 
plasma  conditions  for  which  high-voltage  Orbiter  charging  can  be  expected.  This 
barrier  will  also  suppress  photoemission.  This  appears  likely  to  permit  daylight 
charging  by  this  mechanism,  as  also  was  the  case  for  suppression  by  B  effects. 
We  therefore  predict  that  daylight  auroral-zone  charging  of  wake-region 
dielectric  surfaces  to  potentials  more  than  100V  negative  will  occur  at  least 
occasionally  on  the  Orbiter. 

The  lar  gest  uncertainties  in  our  wake- induced-barrier-effect  predictions  result 
from  three  causes  (Section  4.4).  The  first  is  the  lack  of  available 
measurements  of  electron  and  ion  energy  distributions  made  at  the  same  time 
and  place  with  concentration  measurements  of  individual  ion  species  in  the 
auroral  zone  at  low-orbit  altitudes.  The  second  is  the  difficulty  of  accurately 
predicting  wake-region  eleciron-density  enhancements  caused  by  emitted-electron 
space  charge,  in  the  presence  of  magnetic-field  and  "hopping-conductivity"  effects 
(Sections  3  and  4.3;  see  also  Fig.  3.9).  The  third  cause  is  the  difficulty  of 
accurately  predicting  ion-density  distributions  in  high-voltage  spacecraft  wakes 
(Section  4.3  and  Appendix  C). 


Further  work  is  needed  in  ail  of  these  areas,  and  in  Section  5,  we  have 
presented  preliminary  results  from  a  numerical  simulation  which  was  developed 
in  order  to  make  detailed  investigations  of  the  phenomena  identified  in  Sections  2 
-4.  This  work  involves  calculations  of  floating-potential  distributions  on  infinite 
dielectric  cylinders  in  collisionless  plasma  crossflows  whose  properties  model 
those  of  the  auroral  plasma.  Results  from  this  simulation  indicate  the  presence 
of  multiple  maxima  and  minima  in  the  downstream-side  surface  potential 
distribution,  in  contrast  with  a  monotonic  progression  of  increasingly  negative 
potentials  from  upstream  to  downstream  surfaces,  as  might  have  been  expected. 
These  featur  es  are  caused  by  the  presence  of  ior.  "hot  spots"  or  deposit io-  points 
of  very  beam-like  ion  populations,  which  in  turn  are  the  result  of  a  hig’  ly 
anisotropic  (S^  =  8)  ambient-ion  velocity  distribution  accelerated  into  a  wake 
region  containing  negative  voltages  large  compared  with  the  ion  ram  energy. 

When  combined  with  the  geometrical  complexity  of  a  typical  payload  arrangement 
in  the  Orbiter’s  cargo  bay,  this  promises  to  make  detailed  prediction  of  wake- 
region  surface  potentials  a  formidable  task. 

The  work  presented  in  Sections  2-5  has  involved  initial  explorations  of 
several  distinct  topics  related  to  high-voltage  low-polar-orbit  ehai'ging.  Except 
perhaps  for  our  treatment  of  electron  escape  currents  in  Section  3,  we  have  not 
treated  any  of  tl^ese  topics  exhaustively,  and  our  treatments  therefore  need  in  the 
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future  to  be  extended  in  many  ways.  The  calculation  of  downstream-point  ion 
space  charge  densities  for  a  cylinder,  presented  in  Section  4..'  and  Appendices  C 
and  D,  can  be  easily  adapted  to  spherical  geometry  and  also  to  provide 
downstream-point  ion  current  densities.  These  current  densities  could  then  be 
used  to  provide  simple  estimates  of  required  environmental  conditions  for 
downstream-point  local  charging  in  a  manner  analogous  to  the  estimate  for 
overall  charging  developed  in  Section  2.  Using  available  electron-emission 
yield  data  (Sections  1,3.3,  and  4.5),  these  estimates  should  be  corrected  for 
secondarv-electron  and  backseat tered-elactron  emission  for  use  in  circumstances 
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in  which  it  is  expected  that  these  electrons  can  escape.  The  calculation  of 
emitted -electron  escape  currents  developed  in  Section  3  can  be  readily  extended 
to  provide  corresponding  values  of  surface  currents  due  to  reimpacting  electrons. 
In  fact,  program  SCAPE,  which  is  listed  in  Appendix  A,  already  includes 
provision  for  doing  this,  but  tests  of  this  feature  of  the  program,  conducted  in 
conjunction  with  obtaining  the  escape-current  results  presented  in  Figs.  3.5  and 
3.10-3.17  and  Table  2,  indicated  that  velocity-space  discretizations  which  were 
fine  enough  for  calculating  these  escape  currents  were  not  fine  enough  to  yield 
equivalent  accuracy  in  the  surface-current  calculations.  Surface  currents  need  to 
be  taken  into  account  in  order  to  refine  the  estimates  of  wake  induced-barrier- 
effect  potential  obtained  m  Section  4.  However-,  in  that  situation,  electric  forces 
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on  emitted  electrons  are  directed  toward  rather  than  away  from  spacecraft 
surfaces,  so  the  formulation  of  Sec.  3  is  inapplicable  in  its  present  form.  The 
most  applicable  formulation  would  involve  a  modification  of  the  analytic  surface- 
current  expression  given  by  Eq.  (9)  of  Laframboise  et  al  (1901)  or  Eq.  (6.7)  of 
Laframboise  et  al  (1982b),  to  include  magnetic-field  effects.  The  numerical 
simulation  presented  in  Section  5  should  be  extended  to  include  effects  of 
magnetic  fields,  H^/0+  ion  mixtures,  electron  emission  from  surfaces,  surface 
currents  produced  by  reimpact  of  nonescaping  emitted  electrons,  wake-induced 
barriers,  and  other  spacecraft  shapes,  e.g.  spherical  rather  than  infinite- 
cylindrical,  3nd  effects  of  ion-orbit  curvature  on  space-charge  densities.  All  of 
these  results  need  to  be  compared  with  all  available  experimental  data  on  polar- 
orbit  charging. 

There  remain  two  topics  which  we  mentioned  in  Section  1  but  have  not  treated 
in  this  Report.  The  first  of  these  is  collisional  effects.  We  have  regarded  these 
as  being  inappropriate  for  study  as  part  of  the  work  described  by  this  Report 
because  an  understanding  of  the  important  collisionless  effects  is  a  precondition 
to  understanding  how  these  might  be  modified  by  collisions.  In  any  case,  the 
most  likely  result  of  collisional  effects  on  the  Orbiter  appears  to  be  that 
collisional  ionization  of  fortuitous  gas  releases  may  occasionally  cause 
discharge  of  surfaces  charged  to  large  voltages. 
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The  other  untreated  topic  is  the  possible  modification  of  time-averaged 
spacecraft  wake  quantities  by  nonlinear  effects  of  spontaneous  plasma 
fluctuations.  We  have  briefly  alluded  to  these  effects  in  Sections  3.4  and  5,  and 


we  have  indicated  therein  that  available  evidence  makes  it  seem  unlikely  that 
such  effects  are  important.  However,  both  this  and  the  collisional  problem  are 
largely  unexplored,  and  as  we  have  indicated  earlier  in  this  Section  and 
elsewhere  in  this  Report,  the  various  collisionless  effects  explored  in  Sections 
2-5  are  far  from  completely  studied.  At  present,  it  is  far  from  clear  when 
large  dielectric-covered  spacecraft  such  as  the  Shuttle  Orb  iter  will  begin  flights 
through  the  auroral  zones,  but  the  problem  of  predicting  high-voltage  charging  in 
these  regions  appears  to  be  complex  enough  to  deserve  further  investigation,  and 


this  should  take  place  well  in  advance  of  whenever  such  flights  are  rescheduled 


if  it  is  considered  desirable  to  achieve  detailed  predictive  capability  in  time  for 


them. 
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TABIJE  1.  Low-Earth-Orbit  Conditions** 


!  a..  CHARACTERISTIC  LENGTHS 


i 

j 


j 


Ambient  Debye  Length: 


^  1  cm 


Thickness  of  II  V  and  5kV  sphericaily-symmetric 
sheath***  around  a  sphere  of  radius  3m: 

Thickness  of  ikV  and  5kV  pianar  Child-Langmuir 
sheath: 


2.5m  and  5.8m 


5.0m  and  1 7m 

(WHY  SO  LARGE?) 


Particle  gyroradii: 
ambiem  electrons  (0.1  eV) 
secondary  electrons  (3  eV) 
auroral  electrons  (10  keV) 
ions  (0*;  0.1  eV) 


2  cm 
13  cm 
8  m 


3m  in  "rest"  frame 

27  m*  in  spacecraft  frame 
(larger  inside  a  high-voltage  sheath) 


(b)  CHARACTERISTIC  SPEEDS: 


Ion  thermal  ^oeed: 
Spacecraft  speed: 

Electron  thermal  speed: 


1  km /sec 

^  8  km/sec, 
depending  on  orbit 

100  km/sec, 

for  ambient  electrons 


SITUATION 

IS 

"  MESOTHERMAL" 


*  Particle  motions  do  not  depend  on  the  frame  of  reference  in  which  they  are  viewed. 
The  transformation  from  rest  frame  to  spacecraft  frame  produces  a  V  x  B  electric 
field  ^  0.35V/m  where  V  =  spacecraft  velocity  and  B  =  magnetic  induction;  electric 
fields  inside  a  IkV  sheath  are-'-  200  V/m. 

*v  ceLeeuw  (1967),  p.  1569;  Martin  (1979). 

***  Al’pert  et  al  (1965),  Fig.  72. 
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TABLE  2 


Values  of  the  ratio  1  =  I/I  of  escaping  to  emitted  flux,  for  various  values  of  8, 
the  angle  (in  degrees)  between  the  surface  normal  and  the  magnetic-field  direction,  e, 
the  nondimensional  repelling  electric-fieid  str'ength,  and  the  two  angles  a  and  \p,  also 
in  degrees,  defining  the  electric-field  direction.  The  definitions  of  a  and  \p  are  given 
by  Figure  3.7.  These  four  quantities  appear'  in  Table  2  as  THF.TA,  EPS,  ALPHA,  and 
PSI,  respectively.  For  8  =  0°  and  or  -  0°,  the  electric  field  is  normal  to  the 
surface.  The  results  for  this  case  have  a  numerical  accuracy  of  about  0.2%  or  better; 
see  Sections  3.2  and  3.3.  For  nonzero  values  of  a  and  \p,  the  results  have  a 


numerical  accuracy  of  about  0.4%  or  better;  see  Sections  3.4  and  3.5.  The  results 
given  in  Table  2  are  plotted  in  Figures  3.5  and  3.  10-3. 1  7. 


ESCAPING  FLUXES  AND  SURFACE  CURRENTS 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX 
1  1  11  8  90  4  80 


OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  AUG  29A  .  CASE 

NVY  NVZ 
80  40 


3 


ANCLES  DEFININC  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  0.000  DEGREES ,  PSI 

ESCAP1SC  FLUXES 

EPS  0.00  .20  .50  1.00  2.00  5  00  10  00 

THETA 


=  0.000  DECREES 

20.00 


0.00 

1.000 

1.000 

1.000 

15.00 

964 

.990 

.999 

30  .  00 

.865 

.930 

.977 

45.00 

.706 

796 

.892 

60.00 

499 

.585 

.704 

75 . 00 

.  258 

.311 

.396 

80 .00 

.  173 

.209 

.270 

85.00 

087 

.105 

.137 

37.00 

052 

063 

.082 

89.00 

.016 

.020 

026 

90.00 

0.000 

0  .  000 

0  .000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

.997 

1.000 

1.000 

1.000 

1-000 

.970 

.999 

1.000 

1.000 

1.000 

.856 

.982 

1.000 

1.000 

1.000 

.545 

.802 

.998 

1 . 000 

1 .  000 

.383 

.618 

.968 

1.000 

1.000 

.198 

.341 

.723 

.971 

1.000 

.  119 

.209 

.487 

.810 

991 

039 

.  069 

.172 

338 

.618 

0.000 

0.000 

0 .000 

0.000 

0.000 
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ESCAPINC  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  JUN  13A  ,  CASE  3 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX  NVY  NVZ 

1  1  21  S  90  4  64  64  32 

ANGLES  DEFINING  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  30.000  DECREES,  PSI  =  0.000  DEGREES 

ESCAPING  FLUXES 


EPS 

THETA 

0  .00 

.  20 

.50 

1.00 

2  00 

5.00 

10.00 

20.00 

-60.00 

499 

.527 

.565 

.615 

.676 

.724 

.738 

.745 

-57.00 

.544 

.577 

.623 

.688 

780 

.908 

.981 

1.000 

-54.00 
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.754 
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.818 

.868 
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.976 
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-30.00 

.864 
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.  958 
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1.000 
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2.000 

1.000 

-15.00 

.963 

986 

.997 

1.000 

1.000 
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1.000 

0.00 

1.000 

1.000 

1.000 

1.000 

1 . 000 

1.000 
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1.000 

15.00 

963 
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1.000 
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30.00 
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.932 

.977 

.996 

1.000 

1.000 
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1.000 

45.00 

.706 

.810 

.  9u8 

.976 

999 

1.000 

l  .ooo 

1.000 

60.00 

.499 

615 

759 

.908 

992 

1.000 

1  000 

1.000 

75.00 

258 

.346 

.489 

712 

948 

1.000 

1.000 

1.000 

80 . 00 

.  173 

241 

361 

.579 

892 

1.000 

l.ooo 

1.000 

85.00 

086 

126 

203 

.372 

.730 

.999 

1.000 

1.000 

87.00 

.  051 

.077 

129 

.253 

576 

.992 

1  000 

1.000 

89.00 

.  016 

025 

045 

.099 

278 

1  000 

l.ooo 

90.00 

0.000 

0.000 

0.000 

0.000 

0.000 

0  .  000 

o.ooo 

0.000 

ESCAPING  FLUXES  AND  SJRFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  MAY  30A  ,  CASE  3 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX  NVY  NVZ 

1  1  19  8  90  4  t>4  64  32 

ANGLES  DEFINING  ELECTRIC -FIELD  DIRECTION:  ALPHA  =  60.000  DECREES,  PSI  =  0.000  DECREES 

ESCAPING  FLUXES 


EPS 

THETA 

0  00 

.20 

.50 

1.00 

2  00 

5.00 

10.00 

20.00 

-30.00 

.864 

.874 

.887 

.902 

918 

.929 

.932 

.933 

-29.00 

.873 

883 

.897 

.915 

.937 

.965 

.985 

998 

-28.00 

.881 

,892 

.908 

.928 

953 

.985 

998 

1.000 

-27.00 

.889 

.901 

.918 

.939 

.966 

.994 

1.000 

1.000 

-24.00 

912 

.925 

943 

.966 

.989 

1.000 

1.000 

1.000 

-21.00 

931 

946 

963 

.983 

.997 

1.000 

1.000 

1.000 

-18.00 

.949 

962 

.978 

.992 

1  000 

1.000 

1.000 

1.000 

-15.00 

963 

.976 

.988 

.997 

1  000 

1.000 

1.000 

1.000 

0.00 

1.000 

1.000 

1.000 

1.000 

1  .  000 

1.000 

1.000 

1.000 

15.00 

963 

.983 

.994 

.999 

1 . 000 

1.000 

1.000 

1.000 

30  0" 

864 

.921 

.963 

.989 

999 

1.000 

1.000 

1.000 

45.00 

.706 

799 

.887 

.955 

992 

1.000 

1.000 

1.000 

60 .00 

.499 

.613 

.746 

.879 

.973 

1  000 

1.000 

1.000 

75.00 

258 

.355 

.502 

.702 

.909 

998 

1.000 

1.000 

80  00 

.  173 

.251 

.381 

586 

.848 

.997 

1.000 

1.000 

85.00 

.  086 

.134 

225 

.397 

699 

.990 

1.000 

1.000 

87  00 

051 

083 

.147 

.280 

562 

.971 

1.000 

1.000 

89.00 

016 

.  028 

.054 

.  116 

.284 

.808 

999 

1.000 

90 .00 

0  000 

0.000 

0.000 

0.000 

0 . 000 

0.000 

0.000 

0.000 
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ESCAPING  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  JUL  04A  ,  CASE  3 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX  NVY  NW 

1  1  24  8  90  4  64  64  32 

ANGLES  DEFININC  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  30.000  DECREES,  PSI  =  45.000  DECREES 

ESCAPING  FLUXES 


EPS 

THETA 

0  00 

.20 

-50 

1.00 

2.00 

5.00 

10.00 

20.00 

-67.79 

.377 

.356 

.318 

.242 

.102 

001 

0 . 000 

0.000 

-65.00 

.422 

.409 

.385 

331 

.204 

.011 

0.000 

0.000 

-62.50 

.461 

.457 

.445 

.413 

.321 

.074 

.001 

0.000 

-60.00 

.499 

.502 

.503 

.495 

.451 

.252 

.044 

.000 

-57.00 

544 

.556 

.571 

.589 

.604 

.582 

.500 

.321 

-54.00 

.587 

.607 

634 

674 

.736 

.848 

.945 

.996 

-51.00 

.628 

.655 

.693 

.749 

.837 

.964 

.999 

1.000 

-48.00 

.668 

.701 

.747 

.813 

.907 

.995 

1.000 

1.000 

-45.00 

.706 

.744 

.795 

.865 

.951 

.999 

1.000 

1.000 

-42.00 

.742 

783 

.838 

.906 

.976 

1.000 

1.000 

1.000 

-39.00 

.776 

.820 

.874 

.937 

.989 

1.000 

1.000 

1.000 

-30.00 

.864 

.908 

.952 

.987 

1.000 

1.000 

1.000 

1.000 

-15.00 

.963 

.986 

.997 

1.000 

1.000 

1.000 

1.000 

1.000 

0.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1 . 000 

1.000 

1.000 

15.00 

.963 

.989 

998 

1.000 

1.000 

1.000 

1,000 

1.000 

30.00 

.864 

-93S 

.980 

997 

1.000 

1.000 

1.000 

1.000 

45.00 

.706 

.819 

922 

.984 

.999 

1.000 

1.000 

1.000 

60.00 

.499 

.638 

.803 

.945 

.997 

1.000 

1.000 

1.000 

75.00 

.258 

.397 

.614 

.867 

.992 

1.000 

1.000 

1.000 

80.00 

.173 

.306 

.539 

.834 

.989 

1.000 

1.000 

1.000 

85.00 

.086 

.214 

.471 

.801 

.986 

1.000 

1.000 

1.000 

87.00 

.051 

.182 

.450 

.786 

.985 

1.000 

1.000 

1.000 

89.00 

.  016 

.163 

.430 

.770 

.983 

1.000 

1.000 

1.000 

90.00 

0.000 

.157 

.419 

.761 

.981 

1.000 

1.000 

1.000 

ESCAPING  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  JUL  07A  ,  CASE  3 


NALFA  NPSI  NTH 

NEPS  NTAU 

MINT 

NVX  NVY 

NVZ 

1  1  20 

8  90 

4 

64  64 

32 

ANGLES  DEFININC 

ELECTRIC-FIELD 

DIRECTION : 

ALPHA 

=  60. 

000  DECREES,  PSI 

=  45.' 

ESCAPING  FLUXES 

EPS 

0.00 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

THETA 

-39.23 

.773 

.755 

.724 

.658 

.489 

.074 

.000 

0.000 

-36.00 

.808 

.798 

.781 

.744 

.645 

.279 

.015 

0.000 

-33.00 

.837 

.834 

.828 

.814 

.773 

.587 

.247 

.008 

-30.00 

.864 

.867 

.869 

.872 

.872 

.853 

.803 

.672 

-27.00 

.889 

.896 

.905 

.918 

.938 

.973 

.994 

1.000 

-24.00 

.912 

.921 

.934 

.952 

.976 

.998 

1.000 

1.000 

-21.00 

.931 

.943 

.957 

.975 

.993 

1.000 

1.000 

1.000 

-18.00 

.949 

.960 

.974 

.989 

.998 

1.000 

1.000 

1.000 

-15.00 

.963 

.974 

.986 

.996 

1.000 

1.000 

1.000 

1.000 

0.00 

1.000  1 

.000 

1.000  1 

.000 

1.000 

1.000 

1.000 

1.000 

15.00 

.963 

.984 

.995 

.999 

1.000 

1.000 

1.000 

1.000 

30.00 

.864 

.927 

.971 

.993 

.999 

1.000 

1.000 

1.000 

45.00 

.706 

.816 

.914 

.975 

.996 

1.000 

1.000 

1.000 

60.00 

.499 

.653 

.819 

.943 

.992 

1.000 

1.000 

1.000 

75.00 

.258 

.441 

.690 

.899 

.986 

1.000 

1.000 

1.000 

80.00 

.  173 

.364 

.644 

.882 

.984 

1.000 

1.000 

1.000 

85.00 

.  086 

.292 

.602 

.863 

.982 

1.000 

1.000 

1.000 

87.00 

.051 

.270 

.585 

.854 

.981 

1.000 

1.000 

1.000 

89.00 

.016 

.254 

.566 

.845 

.980 

1.000 

1.000 

1.000 

90.00 

0.000 

.246 

.556 

841 

.979 

1.000 

1.000 

1.000 
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ESCAPING  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 


NALFA  NPSI  NTH  NEPS  NTA'J  MINT  NVX  NVY  NVZ 

1  1  28  8  90  4  64  64  32 


86  JUN  24A  ,  CASE 


3 


ANGLES  DEFINING  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  30.000  DECREES.  PSI 


APING  FLUXES 

EPS 

0.00 

.20 

.50 

1.00 

2 . 00 

5  00 

10.00 

THETA 

-90.00 

.004 

0,000 

0.000 

0.000 

0 . 000 

0.000 

0.000 

-89.00 

.016 

0,000 

0.000 

0.000 

0.000 

0.000 

0.000 

-87.00 

.051 

004 

0.000 

0 . 000 

0  .  OuO 

0.000 

0.000 

-85.00 

.086 

.027 

.002 

0.000 

0 . 000 

0.000 

0.000 

-80.00 

.  173 

118 

.  055 

.009 

0.000 

0.000 

0  000 

-75.00 

.258 

222 

.167 

.089 

.014 

0.000 

0.000 

-72.00 

.308 

.284 

.245 

.178 

.072 

.001 

0.000 

-69.00 

.358 

.346 

.325 

.284 

194 

.028 

0.000 

-66.00 

.406 

.406 

.403 

393 

358 

.222 

.062 

-63.00 

.453 

.464 

.479 

.499 

.529 

.579 

.630 

-60.00 

.499 

.520 

.550 

.596 

.677 

.849 

.973 

-5S.00 

.573 

.607 

.657 

.731 

.848 

.987 

1.000 

-50.00 

.642 

687 

.748 

.832 

.937 

.999 

1.000 

-45.00 

.706 

.758 

.823 

.902 

.978 

1.000 

1.000 

-40.00 

.765 

.819 

.882 

.947 

.993 

1.000 

1.000 

-30.00 

.864 

.914 

.959 

.990 

1.000 

1.000 

1.000 

-15.00 

.963 

986 

.997 

1.000 

1.000 

1  .000 

1.000 

0.00 

1.000 

1.000 

1.000 

1.000 

1 . 000 

1.000 

1.000 

IS  .00 

.963 

988 

.998 

1.000 

1.000 

1.000 

1.000 

30.00 

.864 

932 

.978 

.997 

1.000 

1.000 

1.000 

45.00 

.706 

.809 

.912 

.981 

.999 

1.000 

1.000 

60.00 

.499 

.624 

784 

.936 

.997 

1.000 

l.ooo 

75.00 

258 

.392 

.608 

.871 

.994 

1.000 

1.000 

80.00 

.173 

.311 

.557 

.860 

993 

1.000 

1.000 

85.00 

.086 

.238 

534 

.858 

.992 

1.000 

1.000 

87.00 

.051 

.220 

534 

857 

992 

1.000 

1.000 

89.00 

.016 

.217 

.535 

.857 

.992 

1.000 

1.000 

90.00 

0 . 000 

.217 

.534 

.858 

.992 

1.000 

1.000 

=  90.000  DEGREES 


20.00 


0.000 
0.000 
0.000 
0 . 000 
0.000 
0.000 
0.000 
0.000 
001 
.710 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1  000 
1.000 
1.000 
1.000 


ESCAPING  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  JUN  24A  ,  CASE  6 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX  NVY  NVZ 

1  1  30  8  90  4  64  64  32 

ANCLES  DEFINING  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  60.000  DECREES,  PSI  =  90.000  DECREES 

ESCAPING  FLUXES 


EPS 

THETA 

0.00 

.20 

.50 

1.00 

2 . 00 

5.00 

10.00 

20.00 

-90.00 

.004 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-89.00 

.016 

0  000 

0.000 

0.000 

0  .  000 

0.000 

0.000 

0.000 

-87.00 

.051 

.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-85.00 

.086 

.006 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-80.00 

.173 

.066 

.007 

0.000 

0 . 000 

0.000 

0.000 

0.000 

-75.00 

.258 

.156 

.055 

.004 

0.000 

0.000 

0.000 

0.000 

-70.00 

.341 

.254 

.143 

.037 

.001 

0.000 

0.000 

0 .000 

-65.00 

.422 

.353 

.254 

.122 

.013 

0.000 

0.000 

0.000 

-60.00 

.499 

.449 

.371 

.249 

.079 

.000 

0.000 

0.000 

-55.00 

.573 

.539 

.487 

.397 

.2  27 

.011 

0.000 

0.000 

-50.00 

.642 

.624 

.596 

.544 

.430 

.138 

.004 

0.000 

-45.00 

.706 

.701 

.693 

.676 

.634 

.479 

.228 

.016 

-42.00 

742 

.743 

.745 

.744 

.738 

.699 

.615 

.425 

-39.00 

.776 

.783 

.792 

.804 

.822 

.859 

.901 

.954 

-36.00 

.808 

.819 

634 

.854 

.887 

.947 

.988 

1.000 

-33  00 

.837 

.852 

.871 

.896 

.933 

.985 

.999 

1.000 

-30.00 

864 

.881 

902 

929 

.964 

997 

1.000 

1.000 

-27.00 

.889 

.907 

.929 

.954 

.983 

1.000 

1.000 

1.000 

-15.00 

.963 

.977 

.988 

997 

1.000 

1.000 

1.000 

1.000 

0  .00 

1.000 

1 . 000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

15.00 

.963 

.983 

.994 

.999 

1.000 

1  000 

1.000 

1.000 

30.00 

.864 

.919 

.965 

.992 

.999 

1.000 

1.000 

1.000 

45.00 

.706 

.799 

.896 

.970 

.997 

1.000 

1.000 

1.000 

60.00 

499 

.631 

.793 

937 

993 

1.000 

1.000 

1.000 

75.00 

.258 

.436 

.692 

.914 

990 

1.000 

1.000 

1.000 

80 .00 

.173 

375 

679 

.912 

.989 

1.000 

1.000 

1.000 

85.00 

.086 

.335 

.678 

.910 

.989 

1.000 

1.000 

1.000 

87.00 

.051 

.332 

.678 

.910 

.989 

1.000 

1.000 

1.000 

89.00 

.016 

.333 

.678 

.910 

.989 

1.000 

1.000 

1.000 

90.00 

0.000 

333 

.679 

.910 

.989 

1.000 

1.000 

1.000 

104 


ESCAPING  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  JUL  11A  ,  CASE  3 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX  NVY  NV2 

1  1  21  8  90  4  64  64  32 

ANCLES  DEFINING  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  30.000  DECREES.  PSI  =  135.000  DEGREES 

ESCAPING  FLUXES 


EPS 

THETA 

0.00 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

-90.00 

.  004 

0.000 

0.000 

0.000 

0 . 000 

0.000 

0.000 

0.000 

-89.00 

.  016 

.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-87.00 

.051 

.015 

.001 

0.000 

0.000 

0.000 

0.000 

0.000 

-85.00 

.086 

.  052 

.019 

.002 

0.000 

0.000 

0.000 

0.000 

-82.50 

.130 

.  106 

.074 

.034 

.003 

0.000 

0.000 

0.000 

-80.00 

.173 

.163 

.146 

.115 

.061 

.002 

0.000 

0.000 

-77. SO 

.216 

.219 

.222 

.223 

.214 

.158 

.070 

.006 

-75.00 

.258 

.274 

.297 

.334 

.405 

.598 

.836 

.991 

-72.50 

.300 

.327 

.368 

.438 

.576 

.880 

.997 

1.000 

-70.00 

.341 

.379 

.436 

.531 

.707 

.971 

1.000 

1.000 

-65.00 

.422 

.477 

.558 

.682 

.865 

.998 

1.000 

1.000 

-60.00 

499 

.567 

.661 

.790 

.939 

1.000 

1.000 

1.000 

-45.00 

.706 

785 

.871 

.951 

.933 

1.000 

1.000 

1.000 

-30.00 

.864 

.924 

.969 

.994 

1.000 

1.000 

1.000 

1.000 

-15.00 

.963 

.987 

.997 

1.000 

1.000 

1.000 

1.000 

1.000 

0.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

15.00 

.963 

.987 

.998 

1.000 

1.000 

1.000 

1.000 

1.000 

30.00 

864 

.922 

.970 

.995 

1.000 

1.000 

1.000 

1.000 

45.00 

.706 

.783 

.874 

.959 

.998 

1.000 

1.000 

1.000 

60.00 

.499 

.579 

.691 

.841 

.976 

1.000 

1.000 

1.000 

67.79 

.377 

.454 

.571 

.747 

945 

1.000 

1.000 

1.000 

ESCAPING  FLUXES  AND  SURFACE  CURRENTS  OF  EMITTED  ELECTRONS  IN  UNIFORM  E  AND  B  FIELDS. 

86  JUL  15A  .  CASE  3 

NALFA  NPSI  NTH  NEPS  NTAU  MINT  NVX  NVY  NVZ 

1  1  20  8  90  4  64  64  32 

ANGLES  DEFINING  ELECTRIC-FIELD  DIRECTION:  ALPHA  =  60.000  DECREES.  PSI  =  135.000  DEGREES 

ESCAPING  FLUXES 


EPS 

THETA 

0.00 

.20 

.50 

1.00 

2.00 

5.00 

10.00 

20.00 

-90 .00 

.004 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-89.00 

.016 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-87.00 

.051 

.003 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-85.00 

.  086 

.026 

.002 

0.000 

0.000 

0.000 

0.000 

0.000 

-80.00 

.173 

.126 

.071 

.021 

.001 

0.000 

0.000 

0.000 

-75.00 

.258 

.235 

.202 

.150 

.073 

.003 

0.000 
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Total  Ion  Collection  by  a  Unipotenti* 
Sphere  in  a  Collisionless  Plasma  Flo\ 


Sphere  at  large  negative  potential, 
assuming  that  sheath  thickness 
remains  sphere  radius  r«  , 
using  theory  of  Parrot  et  al  1982/ 
(lower  curve  multiplied  by  X 
1.45)  / 


,l .  _  «i 

7~  'i-t: 


2.69 


1.45 


ion  current 

ion  random 
current 


/  S'  \Sphere  at 

I  /  space  potential, 

/  S'  "any  Debye  length,rT 

XX  ~-8']  ; 

i  =  l7r1//2S:  (ram  current) 

X1  2  1 


0.1  eV  protons 


1  23456789 

an  speed  ratio  Sj  =  U/V 2kTj/mj  - ► 


Figure  2.i.  Dependence  of  ion  current  to  a  sphere  on  Ion  speed  ratio. 
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Figure  2.2.  Dependence  of  spacecraft  surface  potential  on  hot  electron/total  ion 
ambient  flux  ratio.  The  four  circled  points  are  the  largest  values  of  surface  potential 
from  the  numerical  solutions  discussed  in  Sec.  5  and  presented  in  Fig.  5.4. 
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Figure  3.1.  Effect  of  surface  orientation  on  escape  of  emitted  electrons.  In 
(a),  the  spacecraft  surface  is  perpendicular  to  the  magnetic  field  B,  and  the 
emitted  electrons,  which  experience  an  electric  force  -eE  directed  away  from 
tht.  surface,  all  escape.  In  (b),  the  spacecraft  surface  is  nearly  parallel  to  B, 
and  almost  all  of  the  emitted  electrons  return  to  the  surface,  even  though  they 
still  experience  an  electric  force  directed  away  from  it.  Note  that  the 
component  of  E  perpendicular  to  B  results  only  in  an  E  x  B  drift  parallel  to 
the  surface. 
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Ion  shadow  region, 
near -tangential  to  B, 
MOST  LIKELY 
to  charge  to 
large  (negative) 
voltage 


Figure  3.2.  Most  probable  high-voltage  charging  locations  for  a  nonconductive 
spacecraft  in  low-polar-orbit  (auroral-plasma)  conditions. 
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Figure  3.3.  Coordinate  system  for  calculating  electron  escape  fluxes  when  E 
is  perpendicular  to  the  spacecraft  surface.  The  y-coordinate  (not  shown) 
directed  into  the  plane  of  the  Figure. 


igure  3.4.  Example  of  an  electron  orbit  having  zero  initial  velocity.  The 
magnetic  field  B  is  parallel  to  the  (x,z)  plane,  and  makes  an  angle  9  =  75° 
with  the  z  axis,  e  =  i.  Three  gyroperiods  of  the  orbit  {Q  <,  v  <  9n)  are 
shown.  1 1 1 


0  30  ~~ 6a  90° 

e 

Figure  3.5.  Ratio  i  =  I/I0  of  escaping  to  emitted  electron  flux,  as  a 

function  of  the  angle  6  between  the  surface  normal  and  the  magnetic-field 
direction,  for  various  values  of  the  repelling  electric  field  strength  param¬ 
eter  e  =  (E/R)  (Trm/?.kT)^  .  The  electric  field  is  normal  to  the  surface. 
The  result  for  e  =  0  is  given  by  i  -  cos  6.  Values  of  i  plotted  in  this  graph 
are  also  given  in  the  first  section  of  Table  2,  labeled  "A1PHA  =  0.00 
DEGREES,  PSI  =  0.00  DEGREES".  Realistic  values  of  e  for  Shuttle 
high-voltage  charging  conditions  are  in  the  range  30  £  c  £  160  (Sec.  3.2)! 
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Figure  3.7. 


Problem  geometry  when  potential  also  varies  along  surface 
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(b) 


FORCE 

ON 

-eE  =  force  on  SAIL 


Figure  3.8.  (a)  Typical  orbit  of  omitted  electron  when  the  electric  force 

-oF.  on  it.  has  an  outward  normal  component,  but  the  projection  of  -eE 
along  B  has  an  inward  normal  component  (b)  analogous  situation  involv¬ 
ing  sailboat  tacking  into  wind. 
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E 1 ectrone  can  escape. 

Moot  do  if  €  io  1 arge  and 


Figure  3.9.  Dependence  of  electron  escape  and  surface  currents  on 
electric  field  direction  at  surface. 


Figure  3.10.  Ratio  i  =  I/I  of  escaping  to  emitted  electron  flux, 
as  a  function  of  the  angle  6  between  the  surface  normal  and  the 
magnetic-field  direction,  for  various  values  of  the  electric  field  strength 
parameter  e  -  (E/B)  (7rm/2kT)  K  Same  as  Figure  3.5,  except  that  the 
electric  force  vector  -eE  is  no  longer  normal  to  the  surface  (a  is  nonzero) . 
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-90°  -60°  -30 


Figure  3. 12.  Same  as  Figure  3.10,  except  that  -eE  is  no  longer  in  the  same 
plane  as  the  surface  normal  and  the  magnetic  field  vector  is  nonzero). 
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Figure  3. 13.  Same  as  Figure  3.12,  except  that  -eE  is  tilted  further  away 
from  the  surface  normal. 
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-90°  -60°  -30°  0° 


Figure  3.14.  Same  as  Figures  3.10  and  3.12,  except  that  \p  =  90°. 
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90° 


Figure  3.15.  Same  as  Figure  3.14,  except  that  -eE  is  tilted  further  away 
from  the  surface  normal. 
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-90°  -60°  -30°  0°  30°  Q  60°  90° 


Figure  3. 16.  Same  as  Figures  3. 10,  3. 12,  and  3. 14,  except  that  \p  =  135°. 
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-90°  -60°  -30°  0°  30°  q  60°  90 


Figure  3.17.  Same  as  Figure  3.16,  except  that  eE  is  tilted  further  away 
from  the  surface  normal. 
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SPACECRAFT 


SPACECRAFT 


POTENTIAL  MINIMUM 


Figure  4.1.  (a)  General  appearance  of  equipotentials  in  the  near-wake  region 
downstream  of  a  conductive  cylinder  in  a  collisionless  plasma  crossflow  with 
a  large  ion  speed  ratio  S^  and  a  large  ratio  of  cylinder  radius  to  Debye  length. 

(b)  Same  situation  except  that  the  cylinder  is  dielectric  and  emits 
secondary  electrons  and/or  photoelectrons.  In  (b)  the  dotted  curve  is  the 
locus  of  most-negative  potentials  along  surface  normals  (the  negative  "ridge" 
mentioned  in  Sec.  4.2). 
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Figure  4.2.  Normalized  wake  potential  ^  =  -ecfc/kTe^  vs  normalized  distance 
^  =  y/^Qeh  from  the  wake  centerline,  as  given  by  solutions  of  F.qs.  (4.5)  and 
(4.13)  for  various  values  of  the  normalized  wake  half-width  £  given  by  the 
intercepts  of  these  curves  with  the  £  axis.  The  solid  curves  are  for 


normalized  wake  ion  density  n^/n 


ooeh 


equal  to  zero. 


for  nonzero  values  of  this  quantity  as  shown. 


The  dotted  curves  are 


126 


Figure  4.3.  Normalized  wake  potential  well  depth  =  -  e<t>  /kT  u  vs 

°  r  max  max  eh 

normalized  wake  half-width  ^rnax  =  as  given  by  Eq.  (4.9)  for 

normalized  wake  ion  density  nwl  //nooeh  e9ua^  to  zero,  and  by  numerical 
solutions  of  Eq.  (4. 13)  for  nonzero  values  of  this  quantity.  Dotted  curves 


show  limiting  forms  of  Eq.  (4.9)  given  by  Eqs.  (4.10)  and  (4.11). 


127 


I 


Figure  4.4.  General  appearance  of  three  types  of  ion  orbits  impacting  the 
downstream  point  (r  =  r„,  0  =  0)  in  the  presence  of  an  assumed  wake 
potential  distribution  given  by  Eqs.  (4.18)  and  (4.19). 


l/b 


Figure  4.5.  Normalized  ion  density  r).  =  nj/n^  at  the  downstream  surface  point 
(r  =  r  ,0  =  0)  in  the  presence  of  the  potential  distribution  given  by  F.qs. 
(4.18)  and  (4.19),  as  a  function  of  the  ion  speed  ratio  S1  =  U/^kTj/rrn)  , 
for  various  values  of  the  normalized  downstream-point  surface  potential  xo  = 
e<p  /k7\.  Values  of  q^  for  xQ  -  0  are  given  by  Eq.  (4.17).  For  xQ  <  0, 
values  of  q^  are  given  by  Eq.  (Cl  7)  with  T  evaluated  as  described  in  Appendix 
C  and  using  the  computer  program  listed  in  Appendix  D. 
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Figure  5.1.  Expected  general  appearance  of  graphs  of  spacecraft  surface 
potential  as  functions  of  ambient  hot  electron  to  total  ion  flux  ratio  R. 
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EXTRA  GRID  LINES 
AT  85°  AND  95° 


Figure  5.2.  Geometry  of  simulated  spacecraft  used  for  computational 
purposes.  The  outer  boundary  of  the  computational  domain  is  a  square 
whose  sides  are  at  x  =  ±5rg,  y  =  ±5 r  ,  where  rs  is  the  spacecraft 
radius.  The  domain  is  covered  by  the  finite-element  grid  shown  in  Fig. 
5.3.  Within  the  gaps  between  the  conductive  sectors,  the  potential  is 
assumed  to  vary  linearly  with  surface  position  on  the  cylinder.  The  dots 
on  the  conductive  sectors  indicate  locatioas  on  them  where  ion  currents 
are  calculated. 
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Figure  5.3.  (a)  Upper  half  of  finite-element  grid  used  for  computations. 

The  circular  cross-section  of  the  simulated  spacecraft  is  at  the  center  of 
the  grid.  Angular  intervals  in  the  grid  are  10°  except  for  80°  ^  0  ^ 

100°,  where  they  are  5°.  The  radial  intervals  along  the  0°,  90°,  and 
180°  rays  are,  in  order  of  increasing  radial  distance:  0.0004,  0.0016, 
0.0020,  0.0040,  0.0080,  0.0160,  0.0320,  0.1920,  0.7440,  and  3.0. 
The  seven  radial  intervals  nearest  to  the  spacecraft  surface  are  unresolved 
in  (a). 

(b)  Innermost  portion  of  grid  shown  in  (a).  The  first  through 
fifth  radial  intervals  nearest  to  the  spacecraft  surface  are  still  unresolved 


Figure  5.4.  Surface  potentials  as  functions  of  position  for  the  spacecraft 
geometry  shown  in  Fig.  5.2,  for  ratios  R  of  hot  electron  to  total  ion 
ambient  flux  as  shown. 
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APPENDIX  C:  DETAILS  OF  MODEL  CALCULATION  OF 
DOWNSTREAM-POINT  ION  DENSITY  DESCRIBED  IN  SECTION  4.3 

Equation  (4.20)  implies  that: 

v0  ”  ±  mr  ^ C  ”  mq  (Cl) 

We  substitute  this  into  the  equation  of  conservation  of  total  ion  energy  E: 

E  =  ^  mv2  +  q<#>(r,0)  =  \  rnNp2  +  vg2)  +  4  g(©) /r'2  (C2) 

to  obtain: 

r  2  f  c  n* 

V  =  ±  hrr  E  -  (C3) 

r  Lm  L  mr  J  J 

We  introduce  dimensionless  variables  as  follows: 

J  =  L/  (2mkT.)^  r  ;  c  -  C/  mkT.r  2 

IS  IS 

u  =  v(m/2kTi)^  ;  e  =  E/kT.;  X  -  q^/kTj 

g  =  q  g(e)/(kT.rs2)  ;  r  =  r/rs  (C4) 

f  =  (f/nJ(2kT./m) 

00  1 

A  A 

where  f  =  f  (vy,  v^)  is  the  two-dimensional  velocity  distribution  function  given  by: 

f  00 

f  =  d2n/dv  dv  =  f  dv  (C5) 

x  y  z  J 

-oo 


Cl 


Equations  (Cl)  and  (C3)  then  become: 

up  =  ±  (e  -  c/rM  ;  u0  =  ±  [c  -  g(0)]V  r 

The  relation  r  d0/dr  =  v^/v  then  leads  to: 

u  r 


With  g(0)  given  by  Eq.  (4.19),  both  of  these  integrals  can  be  evaluated 
analytically  in  closed  form.  Equations  (4.19),  (C4),  and  (C6)  imply  that: 


(C6) 


(C7) 


g(0)  =  Xo  (1  -  402A2)  <  0  for  101  <  tt/2 


c  =  U0O  +  *0  {C8) 

e  -  u  2  +  y  =  u2  +  c  2:  c 
o  ^o  ro 

2 

From  conservation  of  energy,  we  also  have  e  =  u  ,  where  is  the 
dimensionless  ion  speed  at  infinite  radius,  and  therefore  e  >  0  since  all  ions 

come  from  this  radius.  Since  r  ^  1,  Eq.  (C6a)  then  indicates  that  |ur|  is  an 
increasing  function  of  r,  so  no  turning-points  of  radial  motion  can  occur  outside 

the  spacecraft,  for  those  ions  which  reach  the  downstream  point  (r,0)  =  (1,0). 


However,  in  the  downsteam  region,  g(0)  increases  from  <  0  to  0  as  1 0  J 
increases  from  0  to  tt/2  (Eq.  C8). 


C2 


Therefore,  if  c  <  0,  as  given  by  Eq.  (C6b)  will  become  imaginary  before  1 0 1 
reaches  7r/2.  In  this  case,  there  wili  exist  turning-points  of  circumferential 
motion,  at  angular  positions  9  =  ±  (Fig.  4.4)  given  by  c  -  g(0  )  =  0,  i.e.: 


0  = 
m 


2wue  0/sS-*o 


(C9) 


A  reversed  ion  orbit  having  c  <  0  may  or  may  not  reach  such  a  turning-point 
before  reaching  infinite  radius.  Figure  4.4  shows  typical  ion  orbits  for  c  >  0 
and  c  <  0.  A  numerical  example  of  ion  orbits  having  angular  turning-points,  in 
the  neighborhood  of  a  circular  disk  probe  coplanar  with  a  spacecraft  surface,  has 
been  given  by  Parker  and  Whipple  (1970,  Fig.  5).  Far  from  the  spacecraft,  we 
assume  that  the  ion  velocity  distribution  is  a  drifting  Maxwellian  given  by: 

f  =  (1/tt)  exp  [-  (u  -  S.)2] 

=  (1  /it)  exp  (-ux2  -  uy2  +  2uxSi  -  S.2)  (CIO) 

=  (1/tt)  exp  (-e  +  2uxS.  -  S.2) 

Liouville’s  theorem  for  one  particle  implies  that  f  is  constant  along  each 
collisionless  ion  orbit  (Bernstein  and  Rabinowitz,  1959),  so  we  can  use  (CIO)  to 
find  f  for  any  such  orbit,  specified  by  given  values  of  the  impact  velocity 
components  u^o  and  Uqq  at  the  downstream  axial  surface  point,  if  we  can  find  the 
corresponding  values  of  u  and  u  far  from  the  spacecraft.  For  each  orbit,  (C8) 

X  V 


C3 


can  be  used  to  find  e  and  c.  For  an  ion  orbit  having  c  >  0,  f  can  be  found  by 
first  using  the  integrated  form  of  (C7)  to  find  the  radius  r  at  which  the  orbit 
intersects  the  line  0  -  n/2.  (the  y-axis),  then  using  (C6b)  to  find  u^  =  -Uq,  then 
substituting  this  into  (CIO).  For  c  <  0,  the  procedure  is  more  complicated. 

For  both  c  <  0  and  c  >  0,  we  define: 


G(0)  = 


0 

dO'/lc  -  g(0')]4 

J 

o 


d87(u920  *  4Xq  6'Vir2]4 


for  1 0  I  ^  the  lesser  of  0 


.  -1 
sin 


or  tt/2. 


\2_fTo6 

j  7ru 
1  00 

We  also  define: 


(Cll) 


(  r 


H(r)  = 


'f  /r 


dr  /[  r'(er ' 


For  c  >  0,  we  obtain: 


H(r)  =  (1/VcT ) 


sin  Vc/e  -  sin 


and  for  c  <  0: 


(Cl  2) 


C4 


H(r)  =  -  —  In 


(Ci  3) 


Ve  -  c  +  V-c 
Ve~ -  c/r2  +  \/-c/r 

For  c  <  0,  the  procedure  for  finding  f  is  now  as  follows.  We  first  calculate  r  ^ , 
the  normalized  radius  of  the  first  angular  turning-point  on  the  reversed  ion  orbit 

(Fig. 4. 4).  To  do  this,  we  equate  H(rJ,  from  Eq.  (C13),  to  G(0rn)=7r2/(4V-xo  )> 
from  Eq.  (Cl  1),  and  solve  for  r^.  We  next  divide  H(cc)  -  H(r  J  by  C-(0rrJ  - 
G(-@m)  =  ir2/(2\f-xo)-  The  integral  part  of  this  quotient,  is  the  number  K  of 
complete  continuous  traverses  by  the  ion  orbit  between  -0  and  0m  in  either 
direction.  If  K  is  even,  the  last  angular  turning  point  is  at  0^.  If  K  is  odd,  it  is 

at  In  either  case,  we  next  calculate  bKrj^)  =  H(rp  +  2KG(0m).  If  K  is 

even,  the  "final"  direction  Op  at  infinity,  of  the  reversed  ion  orbit  is  now  given 

by  G(0m)  -  G(0f)  =  H (co)  -  H(r[<+1).  If  K  is  odd,  0p  is  given  by  G(0f)  -  G(-0rr|] 
=  H(oo)  -  Solving  yields: 

ef  =  ±  cos  '-2P-1  (C14) 

[  *'fx^  _ 

for  K  even  or  odd,  respectively,  where: 


(CIS) 


C5 


We  then  have,  at  infinite  radius: 


up  =  -  VF 

ux  =  ur  c°s  0f  (Cl 6) 

u  =  u  sin  &c 
y  r  f 


for  substitution  into  Eq.  (CIO). 

2  2  1 

In  terms  of  the  velocity-space  polar  coordinates  uq  =  (u  +  Uqq  )  ,  a  = 

- 1 

tan  (uqo/,j  ),  anc*  ma^^riS  0 f  hhe  symmetry  of  f  with  respect  to  or,  the  ion 

number  density  at  the  downstream  axial  surface  point  is  now  given  by: 

qo  tt/2 


(Cl  7) 


\  E  ni/nioo  =  2  J  ,  U0dU0  f  !u0’  “•  V  Si>  d“ 


(-x0)2 


=  n.  (y  ,  S.) 
•rV  r 


and,  with  f  evaluated  as  described  in  Eqs.  (Cl)  -  (Cl 6),  this  gives  an  "improved" 

aDDroximation  IcomDared  with  Ea.  (4.17)1  for  values  of  n  .  to  be  substituted  into 
ii  1  1  wi 

Ea.  (4.131  or  (4.14).  The  absence  of  turninp-Doints  of  ion  radial  motion  for  e  > 

j  '  O  l 

0,  which  we  noted  following  Eq.  (C8),  implies  that  all  of  the  ion  orbits  involved 
in  the  integration  of  (Cl 7)  "connect  back  to  infinity"  and  none  of  them  have 
originated  elsewhere  on  the  spacecraft.  Although  a  general  proof  appears 
difficult  to  formulate,  this  probably  means  that  Eq.  (Cl 7)  overestimates  /].  at 


C6 


(r,  6 )  =  (1,0),  whereas  Eq.  (4.17)  underestimates  it.  In  the  applications  treated 
in  Sec.  4.4,  the  spacecraft  geometry  is  closer  to  that  of  Fig.  4.6b  than  Fig. 

4.6a,  and  this  probably  implies  an  even  smaller  ion  density  near  surfaces  in  the 
central  wake  region. 

We  have  evaluated  (Cl 7)  numerically  using  standard  methods.  The  resulting 
values  of  q.  appear  in  Fig.  4.5.  These  values  have  been  computed  to  a  relative 
numerical  accuracy  of  about  0.003  or  better.  The  computer  program  used  for 
doing  this  is  listed  in  Appendix  D. 


C7 


